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We present results for the long-distance asymptotics of correlation functions of mesoscopic one- 
dimensional systems with periodic and open (Dirichlet) boundary conditions, as well as at finite 
temperature in the thermodynamic limit. The results are obtained using Haldane's harmonic- 
fluid approach (also known as "bosonization" ) , and are valid for both bosons and fermions, in 
weakly and strongly interacting regimes. The harmonic-fluid approach and the method to compute 
the correlation functions using conformal transformations are explained in great detail. As an 
application relevant to one-dimensional systems of cold atomic gases, we consider the model of bosons 
interacting with a zero-range potential. The Luttinger-liquid parameters are obtained from the exact 
solution by solving the Bethe-ansatz equations in finite-size systems. The range of applicability of 
the approach is discussed, and the prefactor of the one-body density matrix of bosons is fixed by 
finding an appropriate parametrization of the weak-coupling result. The formula thus obtained is 
shown to be accurate, when compared with recent diffusion Monte Carlo calculations, within less 
than 10%. The experimental implications of these results for Bragg scattering experiments at low 
and high momenta are also discussed. 
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I. INTRODUCTION 

A number of recent experiments have demonstrated the possibility of confining atoms to one dimension. 

This possibility has produced an outburst of theoretical activity over the last years 0, 0, H, EJ fToL ITU Il3l ITiL ITEl 
EU El IB Il9l l2fl I2TL I22I |24| in the field of cold atoms. One important focus of these studies has been on the 
correlation properties of these new one-dimensional systems, with the goal of providing means to characterize them 
experimentally. 

Although one-dimensional models have been the favorite toy of mathematical physicists during the last century, 
they only became experimentally relevant in the late 1960's and 1970's in connection with a number of solid-state 
materials exhibiting very anisotropic magnetic and electronic properties. Later, in the 1980's and 1990's, advances 
in chemical synthesis and nanotechnology made it possible to manufacture materials and devices where electrons 
are mainly confined to move along one or a few conduction channels. The latest experimental developments in the 
field of cold atoms, however, have several distinct features. First of all, the constituent particles are not electrons 
but bosonic atoms (although in the future fermionic atoms can also become available). Second, the "fundamental" 
interaction is no longer of Coulombic type but a short range potential. Furthermore, in the case of dilute atomic 
vapors confined in magnetic and/or optical traps, the degree of controllability over parameters such density and 
interaction strength seems unprecedented. This allows, in principle, to exhaustively explore phase diagrams, or to 
cleanly realize quantum-phase transitions, which so far arc considered as theorists' extreme idealizations of "dirty" 
solid-state phenomena. 

The only a priori limitation offered by the atomic systems, at least from the point of view of a condensed matter 
theorist, seems to be their mesoscopic, rather than macroscopic, size. However, numerical calculations over the past 
decades (e.g. [13 0|) have taught us that many of the behaviors predicted for the thermodynamic limit, already 
manifest themselves at the mesoscopic scale, even when the system under study consists of just a few tens of particles. 
Furthermore, the study of mesoscopic systems is also a brach of modern condensed-matter physics as phenomena 
taking place at the mesoscopic scale have an undeniable interest. Thus the main of motivation of this paper is to 
present a set of tools that can be used to analyze the properties of these mesoscopic one-dimensional (ID) systems. 
Our main theoretical tool in this analysis is the harmonic-fluid approach, which is nothing but the relevant quantum 
hydrodynamics for ID systems. 

The harmonic-fluid approach has a long history pil I21II28. 29], which in some respects culminated with the work of 
Haldane [23, On] . He realized that many one-dimensional models exhibiting gapless excitations with a linear spectrum 
can be described within the same framework. This framework defined a universality class of systems that Haldane 
termed "Luttinger liquids" . The name stems from an analogy with higher dimensional fermionic systems, where the 
equivalent role is played by the (universality class of) Fermi liquids. In the context of one-dimensional Fermi systems 
this approach has a second, more frequently used, name: "bosonization". This refers to the fact that the method 
shows how to describe the low-energy degrees of freedom of the fermions in terms of a bosonic field which obeys a 
relativistic wave equation. However, differently from the Fermi liquids, the class of Luttinger liquids also includes 
one-dimensional interacting boson systems. As we shall discuss below, this has to do with absence of a well-defined 
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concept of statistics in ID. As a consequence, boson systems can display fermion-like properties and vice-versa. One 
well-known example in the field of cold atoms is the behavior of the Tonks gas, where the bosons interact so strongly 
that they effectively behave as free fermions. 

Besides blurring the line that separates bosons from fermions, confinement to one-dimension has another peculiarity 
that is worth discussing. Being all transverse degrees of freedom frozen, fluctuations can only propagate longitudinally. 
This implies that their effect is enormously enhanced. As a consequence, no long-range order that breaks a continuous 
symmetry can exist in the thermodynamic limit, even at zero temperature. 

The harmonic-fluid approach has several advantages over other approaches that are commonly employed to study 
the low-temperature behavior of ID systems. First of all, it is not a mean field theory, and therefore does not break 
any symmetry. Furthermore, it can treat bosons and fermions on equal footing, their difference being manifested 
by the different structure of the correlation functions that describe them. It can also deal with strongly and weakly 
interacting systems at the same time since the low-energy physics parametrized by three phenomenological parameters 
(the particle density and two stiffnesses). These parameters are related to measurable properties of the system, which 
makes the approach conceptually simple. Nonetheless, one has to be aware of its limitations, which essentially are 
related to its "effective field-theory" character. Therefore, the description is a low-energy one, it comes with a built-in 
cut-off, and it is unable to describe the high-energy structure or any model-specific (i.e. non-universal) features. In 
some situations, when the interactions are strong, relating some of the phenomenological parameters to microscopic 
ones can be difficult. But in many such cases one can rely on an exact (i.e. Bethe-ansatz) solution to extract them. 
This will be illustrated here for the case of bosons interacting via a zero-range potential |32l l33j . We will thus be 
able to show that many of the results obtained in the weakly interacting limit from the Bogoliubov- Popov [47L Esj 
approach and its modifications 0,H(j can be recovered and understood within the harmonic-fluid approach. We shall 
also present result for the correlation functions of a number of "toy models" of ID mesoscopic systems, such like a ring 
(i.e. a small system with periodic boundary conditions) and a box (a system with open boundary conditions). Besides 
being analytically tractable, these models are important when comparing with numerical simulations using Monte 
Carlo methods [ ill IT^ . l2ll | . where periodic boundary conditions are used, or the density-matrix renormalization-group 
(DMRG) [34ll35| |. for which open boundary conditions are best suited. Furthermore, in recent times, traps other than 
harmonic have become available 0,H||, which also makes relevant the study of these geometries. 

There already exist some studies of the correlation functions in finite-size systems in the literature, but they have 
mostly focused on fermionic or spin correlations [3^ . l3^.l39ll40l |. The present approach allows to obtain the corrrelation 
functions for both bosons and fermions at the same time, and automatically includes all higher-harmonics, which have 
been omitted in previous treatments. Besides presenting in full detail the generalization of the harmonic-fluid approach 
to the box case, in this paper we also address how to fix the prefactor of the one-body density matrix for bosons 
interacting with a zero-range potential (henceforth referred to as "delta-inteacting bosons"). We find that, when 
properly parametrized, the result obtained by Popov in the weakly interacting limit |48j is accurate even in the st rong 
coupling limit. This is demonstrated by comparing with exact results 01 an d recent diffusion Monte Carlo data plf . 

The organization of this paper is as follows: In the following section, we shall review the harmonic-fluid approach 
following Ref. 30. In this section we also present its generalization to open boundary conditions. In Sect. IIIII 
results for different correlation functions in the box and the ring, as well as finite-temperature expressions in the 
thermodynamic limit, are presented. We consider both bosons and spinless (i.e. spin-polarized) fermions. Sect. IIVI 
specializes the discussion to a system of delta-interacting bosons. This model is relevant for bosonic atoms confined 
to a one-dimensional channel. We show how the parameters needed for the low-energy description can be extracted 
from the exact (Bethe-ansatz) solution, and discuss how to fix the prefactor of the one-body density matrix. Some 
experimental consequences of the correlation functions obtained within the harmonic-fluid approach are presented 
in this section. In particular, the line-shape of the momentum distribution at finite temperature is analyzed as the 
parameters of the system are varied. In Sect.0 we discuss the asymptotic structure of the wave-function of a bosonic 
Luttinger liquid, and obtain an expression for a system with open boundary conditions Finally, the appendices contain 
a proof of a commutation relation used in the main text, an alternative "derivation" of the low-energy Hamiltonian 
using the path integral formalism, as well as a detailed description of how to compute correlation functions using 
conformal field theory methods. Some of the results of this work have been briefly reported in Ref. 0. 



II. THE HARMONIC-FLUID APPROACH 



In this section we are going to review the harmonic-fluid approach in operator language, following the original 
work of Haldane [3(j ■ Some of the results of this section can be also obtained using a coherent-state path integral 
formulation (see Appendix lB")l . 



FIG. 1: Low-energy excitations in a Tonks (= free Fermi) gas. The continuous arrow represents excitations about +pp 
(symmetrically, about — pp) with momenta q « 0. The dashed arrow represents low-energy excitations carrying momenta 
q « ±2pp = ±2npa. 

A. Haldane's construction 

The discussion in this subsection will be independent of boundary conditions, but we shall assume the system to have 
a finite size L. For the most part, the notation is similar to that of Haldane in Ref . l3Ct However, we deviate from it in 
a number of places, sometimes to agree with more recent conventions. We shall work in second quantization most of 
the time. This means that a system of bosons is described using field operators which obey [^(x), w(x')] — S(x — x'), 
and commute otherwise; p(x) = ^(x) 1 $(x) is the density operator. The mean ground state density, po, is fixed either 
by the chemical potential, p, such that po — po(p) (grand canonical ensemble), or by the total particle number in the 
ground state, No, such that po — No/ L (canonical ensemble). Note that the discussion that follows applies to uniform 
systems, the necessary modifications needed to deal with a non-uniform (but slowly varying) ground state density are 
discussed in Sect. Ill El 

As pointed out in the introduction, the effect of long wave-length thermal and quantum fluctuations is enhanced 
by reduced dimensionality. In order to derive a low-temperature description, we need to identify a set of variables 
that describe the low-energy fluctuations of the system. For a bosonic system these variables are the density and the 
phase. At low temperatures, density and phase fluctuations are locally small in the sense to be defined below. To 
give a proper definition of these variables, we need to introduce the phase-density representation of the bosonic field 
operator: 

$ f (i) = V^)e^ (l) - (1) 
Consistently with the bosonic commutation relations ['I'(x), ^(x)] = S(x — x'), phase and density operators obey 

e**^ p{x)e-**^ - p(x) = S(x - x'). (2) 

A convenient way of describing the long wave-length density and phase fluctuations is to split p(x) = p<{x) + p>(x) 
and 4>{x) = 4><{x) + </>>(x), where the p<(x) and </><(x) refer to the "slow" parts (i.e. coarse-grained over distances 
3> Po 1 ) °f the operators, whereas p>(x) and </>>(x) refer to the "fast" or short-wave length parts (see Appendix FBI for 
a more careful definition). In the following we focus on the slow parts, and below, when there is not risk of confusion, 
we will denote as <j>(x) the slow part of the phase operator (i.e. 0<(x)). Furthermore, since the density fluctuates 
at low temperatures about the ground state value, po, it is convenient to introduce the operator II (x), defined by 
P<(%) = Pa + n (x). 

We shall first consider the phase fluctuations. In Appendix IaI we show that the slow parts 4>{x) and II(x) are 
canonically conjugated fields, i.e. 

[II(x),0(x')] = id(x-x'). (3) 

It is worth pointing out that this commutation relation holds, for the slow parts of the density and phase, in arbitrary 
dimensions. However, since in ID fluctuations are constrained to propagate on a line, this causes the aforementioned 
enhancement of their effect. This makes it impossible, in the thermodynamic limit, to define of an order parameter. 
Thus we expect that (e 4 "^)) = in ID, in contrast with the situation for d > 1, where this operator can acquire a 
non-zero expectation value (i.e. there can be off-diagonal long-range order). Indeed, in one dimension (e 1 ^^) = 
also in finite-size systems (see appendices [C] and [DJ) . 

We next consider long wave-length fluctuations of the density, which have been parametrized as po + n(x). n(x) 
describes locally small fluctuations of wave-length 3> p^ 1 . However, this does not suffice to describe all possible 
low-energy density fluctuations. It is thus necessary to distinguish between "low-energy" and long "wave-length" ; 
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There can be low-energy density fluctuations with wave-length ~ p$ 1 or shorter. To see this in a non-trivial case, 
it suffices to consider the case of the impenetrable-boson (Tonks) gas. It is known [33l that the excitations of 
such a system are those of a free fermi gas of the same density. In a Fermi gas the momentum of the fastest particle 
(the Fermi momentum) is related to the density by the formula pp = ttpq. There are two Fermi points (see Fig.^l, 
corresponding to p — ±pf- The low-energy and long wave-length density fluctuations of this system are the small 
momentum particle-hole excitations around each one of the Fermi points. However, a fermion can also be excited 
with very low energy from one Fermi point to the other, thus producing a density fluctuation that oscillates as 
cos2ppx = cos27rpox. The long wave-length density fluctuations described by n(a;) lead to small changes in the local 
Fermi momentum: pp{x) = pf + ttII(x). This affects shorter wave-length density fluctuations, which now oscillate as 
cos 2 f x cIx'pf(x') = cos20(x), where we have introduced the auxiliary field <d(x), related to II(x) by means of the 
following expression pl.l30|: 

-d x Q{x) = p +U(x). (4) 

Note that the integral on the left hand-side of Q from x = to x = L must be equal to total particle number operator 
N, namely 

Q(L) - 6(0) = ttN. (5) 

This is an important topological property of 0(x). It tells us that on a global scale the changes in Q(x) are related 
to changes in the total particle number. Hence, the configurations of Q(x) can be regarded as functions that increase 
monotonically from x — to x — L (see Fig. [5J ■ Thus it is tempting to associate the location of particles with the 
points where &(x) equals an integer multiple of tt. Physically this amounts to trading discrete particles by solitons 
or kinks in Q(x) (see Fig. |2J. The construction of a low-energy projection of the full density operator, p{x), which 
reflects the discrete nature of the particles and therefore can describe the shorter wave-length low-energy fluctuations 
follows from this idea [3(| . First one sets: 

+oo 

p(x) = d x e(x) J2 S(Q(x)-mr), (6) 

n— — oo 

which in the particle-soliton sense is equivalent to the first quantized form of the density operator, p(x) — yV —1 S(x— Xi) 
provided that Xi is interpreted as the position operator of the particle-soliton and one uses <5[/(a:)] = 5(x — xo)/\f'{xo)\, 
where f(xo) = 0. The previous expression can be rewritten in a more useful way with the help of Poisson's summation 
formula: 

+oo +oo <.+oo 

£ f(n)= Yl / dzf{z)e^\ (7) 



A 0( X ) 




FIG. 2: Typical configuration of Q(x). Its value changes by it every time x surpasses the location of a particle. Note that Q(x) 
is a slowly varying field. 
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which yields: 

+00 +00 



p(x) = -^6(x) V e 2 "" 6 ^ = [p + H(x)] y e 2mie ^\ (8) 



m— — 00 m— — 00 



This is the sought representation of the density operator. The m = term is precisely po + n(x), and describes 
the long wave-length fluctuations (i.e. those with momenta \q\ <C po), the m = ±1 terms describe fluctuations with 
q w ±2npo, while m = ±2 those with q w ±4-7rpo, etc. (see Fig. EJ). 

We next take up the construction of the boson field operator ^(x) in terms of <d(x) and 4>(x). According to Eq. QJ, 
this requires finding a representation for the square root of the density operator, Eq.©. At this point it is useful to 
recall Fermi's trick: [£(#)] 2 = A 5(6), where the constant A depends on the particular way the Dirac delta function is 
defined. Extracting the square root yields y/I{0) = A" 1 / 2 8(0). Thus, using Poisson's formula Q and multiplying by 
e -uf>{x) f rom the right, one arrives at 30]: 

+00 

¥(x) ~ [ Po + n(x)}^ e 2mie[x) e' 1 ^ . (9) 

m — — oo 

The symbol ~ means that the field operator is given by the expression on the right up to a prefactor. This prefactor 
is not determined independently of the way we choose to exclude the high-energy fluctuations (i.e. the fast modes) 
from the low-temperature description. This is usually done using a cut-off in real or momentum space, and different 
cut-off schemes lead to different prefactors. 



j© (x) 




FIG. 3: Effect on the string operator e 1 ®^ of the creation of a particle at two different points x\ (a) and X2 (b). When a 
particle is created (or removed) at Xi, the string operator at X2 < xi is not affected (a). However, when the particle is created 
(or removed) at X2 the value of the string at xi picks up an extra factor of — 1. Since the operators ^ F (x) contain only odd 
powers of e lS x ' , they are sensitive to these changes, which implies that they anti-commute. 



The construction of a fermionic (i.e. anti-commuting) field operator is also possible. One only needs to realize that, 
since @(x) jumps by tt every time a particle is surpassed (Fig. [3J) , the string operator e ,e ^' alternates between +1 
and —I (see Fig- EJ - This implies that 

+00 

*t,(a:) =tft(a.) e <eO») „ [ Po + n(x)}^ ^ e (am+i)ie(») e -i<f>( x ) ( 10 ) 

m— — 00 

anti-commutes at different positions. To see this, consider the product 1 >S? F (x2)'& F (xi) assuming that xi > x%, for 
instance. Since the operator ^ F (xi) acts first by creating a particle at x\, it does not have any effect on the string 
operator at x-i < x\ (see Fig. |3 case (a)). On the other hand, if we instead consider ^ F (xi)^/ F (x2) (%i > x 2 ), the 
operator ^ F (x2) acts first, thus creating at x-i < x\. Therefore, the string e ie< - Xl ^ of the second field operator picks 
up an extra minus sign (see Fig. case (b)) relative to the case where a particle is created at x\ first. Thus the 
operator defined by i|IU|) anti-commutes at different locations, and therefore describes fermions instead of bosons. The 



6 



same conclusion can be reached for a product of two annihilation (or one creation and one annihilation) operators at 
different points [73 ]. 

The explicit construction of the low-energy representations of commuting (i.e. bosonic) and anti-commuting (i.e. 
fermionic) fields shows how this approach treats bosons and fermions on equal footing. On physical grounds, the fact 
that in ID transforming bosons into fermions is possible is not a surprise. The reason |30| is that when one tries to 
exchange two interacting particles (or elementary excitations) in ID they must necessarily collide, and therefore the 
statistical phase cannot be separated from the phase shift associated with the collision (see Fig.0J. 



B. Low-energy effective Hamiltonian and momentum operators 



So far we have considered the kinematics of the low-energy description of a ID quantum fluid. In other words, 
we have introduced the variables that describe low-temperature states of the system and found their relationship to 
density and field operators. The next step is to consider the dynamics, that is, the Hamiltonian. The starting point 
for our considerations is the following Hamiltonian for bosons interacting via a general two-body potential, v(x) |73j : 



h 2 f L 1 f L 

H = — j dxd x ^{x)d x ^{x) + -J dx dx v{x - x')p{x)p{x) 



(11) 



Before proceeding any further, a number of comments arc in order: The above Hamiltonian is assumed to describe the 
situation where all particles lie in the lowest level of a transverse confining potential. This happens when the chemical 
potential, fx, is smalller than the transverse confinement energy, huj±. Regarding the longitudinal confinement, we will 
assume in the following that it is either absent (case of periodic boundary conditions) or that it can be approximated 
by two infinite barriers placed at x = and x — L (case of open boundary conditions) . The effect of a smooth potential 
in the longitudinal direction will be discussed in Sect. Ill El By restricting ourselves to the lowest transverse level, 
we assume the system to be effectively one-dimensional. If Uq (y, z) denotes the lowest transverse orbital, the three 
dimensional boson field operator can be written as ^^(r) = ^(x) u^(y,z) + ^(r), where ^(r) describes bosons 
in the higher energy transverse levels. At low temperatures (i.e. T < fi < Tilo±), only virtual transitions to higher 
transverse levels are permitted, which lead to a renormalization of the interaction potential v(x) (see e.g. Ref. |fj) 

To obtain the low-energy effective Hamiltonian we use the operator identities derived in the previous section keeping 
only the leading terms, which are quadratic in the gradients of the slowly varying fields <d(x) and 4>(x) (alternatively, 
one can linearize the equations of motion for the density and the phase-gradient in terms of the gradients <fi and 0). 
The result can be generally written as 



vj (d x (j)(x)) + v N (d x Q{x) - irpo) 



(12) 



The history of this Hamiltonian goes back to the pioneering work of Tomonaga |26| on one-dimensional electron gases. 
He wrote it in a very different way, but as we shall see in the following sections, both forms are essentially equivalent 



time 





space 



FIG. 4: When two particles are exchanged in one-dimension to test their statistics, they cannot avoid collision. Thus the 
statistical phase and the scattering phase shift cannot be separated. 
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and describe the same collection of harmonic oscillators, whose quanta, the "phonons" correspond to low-energy 
density and phase fluctuations. The first application to a bosonic system was done by Efetov and Larkin psj . who 
considered a system of tightly bound electron pairs (the "BEC" limit of a ID superconductor) and wrote H g in a 
form similar to Eq. 1)120. The form used here is due to Haldane |30tl3l| who, building upon and extending the work of 
Tomonaga j2||, Lieb and Mattis [2?|], Luther (2^, Efetov and Larkin [2^,... introduced the concept of (Tomonaga-) 
Luttinger liquid as a universality class of ID systems with gapless, linearly-dispersing, excitations. 
For an interaction whose range R^> po , one should replace the second term in H e s by 

^2 J dx dx' v(x - x')(d x <d(x) - np )(d xl <d(x') - irp Q ). (13) 

However, provided that the Fourier transform of the interaction potential v(q) — J Q L dx e~ lqx v(x) is not singular as 
q — ¥ 0, Ea. H12|) will hold at sufficiently low energies, and below we shall assume that this is indeed the case (a notable 
exception is the Coulomb interaction). The term proportional to (d x Q — Trpo) 2 in Eq. i|12|) thus corresponds to the 
interaction energy of the long wave-length density fluctuations. The density stiffness vn has dimensions of velocity. 
In the following subsection, it is shown that it is inversely proportional to the compressibility of the fluid. The phase 
stiffness Vj has also velocity units and it is found to be proportional to the superfluid fraction. Therefore, both 
parameters must be regarded as phenomenological. They can be extracted from an exact solution of the microscopic 
model (when available), from numerical calculations or, ultimately, from experimental data, which still allows one to 
correlate the results from different experiments. 

For many models on the continuum (i.e. those which do not require a lattice to be defined) the term proportional 
to (d x (j)) 2 in Eq. (|12|) follows from the kinetic energy operator after using Eq. @, i.e. 

I! ' < L dxd x ^(x)d x ^(x)^^- [ L dx {d x cj>{x)f. (14) 



2M J * v ' * v ' 2M ... 

Hence Vj = hnpo/m = vp, where vp, is the Fermi velocity of a gas of free spinless fermions of density po- As we show 
in the following subsection, this relationship between vj and vp is ensured by Galilean invariance, and will not hold 
if this symmetry is broken. This is generally the case of lattice models, such like a ID system of bosons hopping in a 
sufficiently deep optical lattice 

Another interesting operator is the total momentum, which can also be expressed in terms of Q(x) and 4>{x). The 
derivation is similar to that of the Hamiltonian, starting from the second quantized form, 

P =%\ L [*Hx)d x *(x) - d x &(x)9(xj\ , (15) 
^ Jo 

we keep the leading terms in the gradients of 0(x) and <fi(x), and obtain 

P = - [ d x e(x)d x ct>(x). (16) 



o 



We close this subsection with a remark about notation. It is customary in the literature on bosonization to work 
with the field 

8(x) = 6(x) - np x, (17) 



instead of O(x). Hence, H(x) = d x 6(x)/ir. It is also common to introduce the parameters K — \/vj/vn and 
v s = y/VNVj, such that 



H cS = ^ / dx 
* Jo 



- {d x ^{x)) 2 + I (H(x)) 2 

7T K 



2tt Jo 



dx 



K (d x (b(x)) 2 + 1 (d x 8(x)f 



(18) 



As we shall see shortly, v s is the phase velocity of the low-energy excitations (sound waves). However, the dimensionless 
parameter K is related to the strength of quantum fluctuations (see Appendix |B|) . 



C. Particles in a ring: periodic boundary conditions. 



Our first task in this subsection will be to find appropriate mode expansions for the fields Q(x) and 4>[x), such 
that the Hamiltonian, Eq. i|12[l . is diagonalized. Before doing it, we need to find out how to implement the boundary 
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conditions in terms of these fields. Therefore, consider a system obeying periodic boundary conditions. As commonly 
introduced in the literature, this seems a mere mathematical convenience that simplifies the calculations before taking 
the thermodynamic limit. However, nowadays there exist experimental realizations of these BC's. In the case of cold 
atoms one can think of a quantum degenerate atomic gas in a tight toroidal trap [3(|. Thus, if the boson field obeys 
+ L) = *t(x), then equations and (1SI9P imply that: 

<d(x + L) = Q{x)+7rN, (19) 
(j)(x + L) = (f>{x)+nJ, (20) 

where N is the particle-number operator, and J is an operator whose eigenvalues are even integers so that (— \) J = +1. 
These are very important topological properties of the phase and density fields. In particular, they imply that the 
system has states that can be labeled by the eigenvalues of J and TV. We show below that N and J are conserved 
quantities (i.e. they commute with H e g) and can be used to label topologically excited states of the system. In the 
case of N this is not surprising because the microscopic Hamiltonian, Eq. iJHJ, conserves the total particle number. 
However, we will see that J is associated with the possibility of quantized persistent currents. 

We next write down the mode expansions for Q(x) and (f>(x), which obey and l(2TJ|) . as well as the commutation 
relation: 



[d x Q(x),<f>(x')] = inS(x - x'). (21) 



The appropriate expressions read: 



Q(x) = o + ^iV+i^|^|%-l^ 2 [ e ^6( (Z ) + e-^t ((z) ] j (22) 

q^O q 

<t>{x) = o + ^J + i^|^|%^M/ 2 sgn( g )[e^6( g ) + e -^6t(g)]. (23) 
q^a q 

The operators b(q) and b^(q) have the commutation relation: 

[b{q),tf{q')]=6 q<q ,, (24) 

commuting otherwise. Momenta are quantized as q — 2irm/L, where m — 0, ±1,±2,±3, . . . The pairs (N,cj)o) and 
(J, 9q) are conjugate action-angle variables which obey: 

[N^e-^ ] = e -i *°, (25) 
[J,e- l6 °] = e' l6 ° (26) 
[J,N] - [4> ,9 ]=Q. (27) 

Introducing the expressions l|22|) and (|23|l into Eq. (|12fl . one obtains: 

H cS = M<7) bHq)b(q) + g|(7V - iV ) 2 + J 2 + const., (28) 

where cu(q) = v s \q\ for q <C p^ 1 ■ It now becomes clear that this Hamiltonian describes collective phonon-like excitations 
(sometimes called Tomonaga bosons 26]), which disperse linearly in the long wave-length limit |g| — > 0. This result 
is perhaps not striking if one deals with bosons because it is already obtained from Bogoliubov's theory (althou gh in 
ID this theory is inconsistent in several respects, which can be cured only in the weakly interacting limit |sl I20ll4^| ) . 
However, for interacting fermions it has striking consequences, because it means the absence of individual (i.e. particle- 
like) excitations in the low-energy spectrum, and the break-down of Fermi-liquid theory in ID. 
The expansions (|22|l and ill'MI) diagonalize the momentum operator as well: 

P = ^NJ + J2^(q)b(q). (29) 

q^O 

As O(x), <j){x), H c ff and P describe long wave-length fluctuations, the momentum sums must be cut-off at a momentum 
a -1 < min {-R -1 , (? c }. For short range interactions, a -1 < q c ; q c is then fixed by demanding that hv s q c = /z, i.e. it is 
an estimate of the momentum where the excitation spectrum deviates from the linear behavior. 
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FIG. 5: Typical low-energy excitation spectrum of a Luttinger liquid with periodic boundary conditions. The shaded area 
corresponds to the continuum of multi-particle-hole excitations. The Hamiltonian H e ff, Eq. 12811 describes the excitation 
spectrum for q « nJpo, with J = 0, ±2, ±4, . . .. Near q m the excitations are linearly dispersing phonons: cu(q) ~ 7n> s |5| (i.e. 
superpositions of particle-hole excitations about the same Fermi point, see Fig. 0. For q ~ nJpo, where J — ±2, ±4, . . ., the 
excitations involve the creation of J current quanta in the ring. For instance, J = ±2 correspond to the low-energy transitions 
between Fermi points shown in Fig. [5] (dashed arrow). For delta- interacting bosons, this picture holds in all regimes, and for 
the particle- hole spectrum agrees with the spectrum obtained from Bethe ansatz in Ref. |5ll |. 



From the previous expression for the Hamiltonian the (inverse) adiabatic compressibility at zero temperature can 
be obtained: 



-1 _ 2 



pIl 



d 2 E (N) 



dN 2 



N=N n 



hnv s pl 



(30) 



where E (N) = (H)n is expectation value of the Hamiltonian, Eq. I|28[l . taken over the ground state with N particles; 
hence u/v oc Kg 1 . Using the previous expressions, we can also write: 



1 



v N 







\d 2 E Q (Ny 








dN 2 


N=N 



(31) 



which will be useful in extracting the density stiffness from the Bethe-ansatz solution of the delta-interacting bosons 
in Sect. EH 

We next find the relationship between the phase stiffness, vj, and the superfluid fraction, which we denote as p s 
below. To this purpose, we consider twisted boundary conditions instead of PBC's: 



VHx + L) 



- fa *t( 3 



(32) 



Therefore, the phase field obeys the modified boundary conditions: </>(x + L) = <f>{x) + ttJ + a, such that (— 1) J = +1. 
This means that in Eq. we have to shift J — > J + a/ir. As a result, the ground state energy and momentum are 
also shifted (recall that in the ground state (N) = iV and (J) = 0): 



(-Heff)a^O — (H c ff) a= o 

(P) a ^o - (P) a =o 



i fhv 



2 

TIN, 



L, 



(33) 
(34) 



Thus we see that the system responds to the twist in the BC's by drifting as whole with constant velocity v s t = 
((P)a^o — (P) a=o) I {M No) = ha/ ML. The superfluid fraction, p s , can be obtained by regarding the shift in the 
ground state energy as the kinetic energy of the superfluid mass M{p s L): 



By comparing (|33|) and l|35|) we can identify 



ft Ps /a\ 2 



(z)' 



L. 



hvj 



n ps 

M 



(35) 



(36) 
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Otherwise, in view of Eq. (|33[1 . we could have defined the phase stiffness in an analogous way to the density stiffness: 

ttL f d 2 E (a) \ 

vj= ^{-a*-)^> (37) 

where E (a) = (H) a ^ . Hence, vj is seen to be related to the response of the system to a phase twist just like vm is 
related to the response to a change in particle number. For charged particles the phase twist can be thought as the 
result of a magnetic flux that threads the ring [57j . 

We next specialize to Galilean invariant systems. As we have already anticipated in the previous section, for these 
systems vj = vp = hnpo/M, and therefore p s = po- To prove this identity we perform a Galilean boost where the 
particle momenta are shifted: pi — > p[ = pi + TiQ. The Hamiltonian and momentum operators transform as: 

P' = P + hNQ. (39) 

Therefore, the ground state energy and momentum are shifted: E' — E = N(TiQ) 2 /2M, and Pq = KNQ. Upon 
comparing the latter expression for P with Eq. i|29|) . we see that, to leading other, Q = irJ/L. Hence, E' Q = 
Eq + (h7r) 2 p J 2 /(2L). However, from Eq. J2HJ, the energy of the boosted state is E' = E + Tittvj J 2 /(2L). From 
these last two expressions it follows that vj = vp. 

The previous analysis also allows to understand why no term involving the product d x @(%)d x <fi(x) appears in the 
low-energy Hamiltonian, Eq. 1)12(1 . The reason is that this product is the momentum density (cf. Eq. pfifl). Therefore, 
the presence of such a term in the Hamiltonian would indicate that we have not chosen the reference frame where the 
system is at rest. Thus one can get rid of it by a suitable Galilean transformation (cf. Eq. (J3HJ). 

We end this subsection with a discussion of the selection rules for the eigenvalues of J and N. For bosons with 
PBC's we have found the following selection rule: (— l)' 7 = +1, i.e. J has eigenvalues that are even integers. However, 
obtaining the selection rule for fermions becomes somewhat messy using the above methods. To derive the selection 
rules for both bosons and fermions we recall here a construction due to Mironov and Zabrodin [Hq| . The construction 
works directly with the many-particle wave function and therefore does not rely on the previous formalism. Consider 
a rigid translation of the system by a distance a, i.e. 

e lPa/h <f>(x 1 ,x 2 , ...,x N ) = + a, x 2 , ■ ■ ■ , x N + a) (40) 

Let $ be an eigenstate of the momentum operator P with eigenvalue P and set a = L/N = Pq 1 , x\ = x,x 2 = 
x + a, . . . , xn — x + (N — l)a. Thus, 

e lPa/h ^(x, x + a, . . . , x + (N - l)a) = $(x + a, x + 2a, . . . , x) = (il)^ 1 *^, x + a, . . . ,x + (N - l)a), (41) 

where the plus sign corresponds to bosons and the minus to fermions. The last equation follows from permuting the 
particle coordinates in the translated wave function to sort them in increasing order. The operation involves N — 1 
transpositions and therefore for fermions the signature of the permutation is (— l)^ -1 . If we choose the wave function 
to be an eigenstate of P with P = hirJN/L, we arrive at the selection rule: 

i-iy = (±i) N -\ (42) 

which reduces to (— 1) J = +1 for bosons, and to: 

i-iy = -(-i) N (43) 

for fermions. In the thermodynamic limit these selection rules are not very important because one can effectively 
treat N and J as having continuous eigenvalues. Furthermore, the terms in the Hamiltonian 128|) proportional to 
(N — Nq) 2 and J 2 are of order L _1 and disappear as L — > oo. However, for the mesoscopic systems which concern us 
here, the selection rules can be important in certain situations where one needs to keep track of finite-size effects and 
the discreteness of particles. 



D. Particles in a box: open (Dirichlet) boundary conditions. 



We now turn our attention to systems with open (Dirichlet) boundary conditions (OBC's). In cold atom systems 
these conditions have been already experimentally realized (in an approximate way) |J] using a microchip trap where 
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the potential was shaped to a square-well with very high barriers, which can be approximated by perfectly reflecting 
walls. Less restrictively, one can assume that there are two points, x — and x — L, where current density vanishes, 
and this corresponds to the equilibrium state in the experiment of Ref.|3, where 4 He was confined in a long nanopore. 
To implement the latter BC's, we first note that, from the continuity equation, 

d tP (x,t) + d x j{x,t)=0, (44) 

and using that p(x,t) ss d x O(x)/ir, it follows that the current density j{x, t) w — dtQ(x, i)/7r. Hence, demanding that 
j(x = 0) = amounts to dtQ(x = 0, t) = 0, i.e. <d(x = 0, t) = Ob = const. However, what is often understood by open 
boundary conditions is something more restrictive than this: it is demanded that p(x = 0, t) = 0. This is achieved 
by the further requirement that 9 b ^ 0, ±7r, ±2-7r, ... In other words, Q(x — 0) must be pinned at real number that is 
not a multiple of ir. To understand this, we need to go back to Eq. JfjJ), 

p(x) = d x Q{x) ^2 S(e{x)-mr). (45) 



From this expression, one can see that p(x — 0) — provided that Q(x — 0) = 8b 7^ nir, where n is an integer. If 
p(x = 0) vanishes so does W(x — 0) cx -J 'p(x = 0). What about the other end of a finite system (i.e. x = L)l The 
property 

[Q(L) - 6(0)] = Ntt (46) 

also fixes &(L) = 0(0) + Nir = 9b + Ntt. Thus the boundary condition is automatically satisfied at x = L. 

We now find the appropriate mode expansions for 0(x) and 4>{x). The requirements are the same as in the previous 
subsection, namely that the BC's and the commutation relation of d x Q(x) and 4>{x) 1 Eq. I|21(l . must be fulfilled for 
< x,x' < L. Thus we arrive at the following expressions (remember in what follows that 9b is just a real number, 
not an operator): 

e(x) = 9 B + ^N + iY,(^f e- a «' 2 sm{qx)[b{q)-b\q)l (47) 

q>0 q 

<t>(x) - ^0 + ^ (-^) ' e-^' 2 cos(qx) [6(g) + b\q)] (48) 

where a -1 < min{i? _1 ,g c = p/hv s } and q — rmr/L, with m = 1,2,3, .. . Notice that in this case only one pair of 
action-angle operators is needed, (N,(j)o)- This is because in a box, as opposed to a ring, there cannot be persistent 
currents, and therefore J is not a good quantum number. 

The above mode expansions diagonalize the Hamiltonian and render it as follows: 

H cS = Ml) bHq)b(q) + ^(N - N f , (49) 

q>0 

with u>{q) = v s q for q <^ pq. The restriction q > means that q cannot be interpreted as the momentum of the 
excitation but as its wave number. The sound waves in a system with OBC's are standing waves and therefore do 
not carry momentum. Consistently, the momentum operator vanishes, i.e. P = [Fl ]. 



E. Effect of a slowly varying confining potential 



In most of current experimental setups cold atoms are confined in harmonic traps. In this section we are going to 
discuss how the harmonic-fluid approach must be modified when a smooth potential is applied in the longitudinal 
direction. We shall distinguish two situations. In the first we consider that the external potential is weak and show 
that much of what has been said above applies with small modifications. However, the harmonic potential does not 
belong to this class. In other words, it is not weak, but this does not mean that the harmonic-fluid approach cannot 
be adapted to this case. However, calculations are no longer analytically feasible (except in certain limits). Some 
results for harmonically trapped gases of bosons |l4j ] and fermions [l5| are already available in the literature, but 
since in general no explicit expressions for the correlation functions are available it is hard to extract much qualitative 
information from these results. In the Tonks limit, one can rely on the fermion-boson corresondence established by 
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Girardeau |49( , which also allows to make a beautiful connection with the theory of random matrices in the gaussian 
unitary ensemble (GUE). This connection makes it possible to obtain the density profile [41j, the form of the Friedel 
oscillations |H3] as well as asymptotic forms for the one-body density matrix 0] . It is also worth mentioning that for 



fermionic systems one can use constructive bosonization |31j . and some analytical results can be thus obtained 22] . 
However, the problem with this approach is that the interactions between cold fermions do not have a simple form in 
the harmonic-oscillator basis. Therefore, it becomes hard to assess the effect of the different matrix elements on the 
low-energy physics. 

Introducing an external potential amounts to adding to the Hamiltonian in Eq. Qllfl the following term 



dx u ext (x)p(x) 



(50) 



Provided that u cx t(x) is weak (i.e. |it ex t(aO| *C p) and varies slowly over distances on the scale of a, it couples only to 
the slow part of the density operator, namely to d x Q(x)/-K. Hence, 



-Hcff 



dx 



~2~7 



K {d x <j>{x)) + -1 (d x Q{x) - vrp ) 2 

K 



Wext(z) 



d x @(x) 



Up to a constant, H e g can be rewritten as: 

r-L 



~2~7 



dx 



n 



K{d x( t>{x)y 



K 



where pa(x) — po + ^Po(x), and 



Sp (x) 



K 



{d x <d(x) - irp (x)Y 



dp\ 



Ucxt(x) = - ( —j u ext (x). 



(51) 



(52) 



(53) 



The last expression follows from the relationship found in Eq. (|3U|I between K/hnv s and the compressibility. The 
result from linear response theory is thus recovered, since M ex t(^) enters the Hamiltonian as (minus) the chemical 
potential. In order to render the Hamiltonian to its form in the absence of u ex t it is convenient to shift 

(54) 



(55) 



G(x) =Q{x)-tt I dx' [p (x')-p ], 
such that the second term in Eq. I|52|l becomes: 



2nK 



dx 



(d x e(a 



The expansions in modes given in previous sections can now be used to diagonalizc H c g with &(x) playing the role of 
0(x). This leads to the same form for the spectrum, but the shift (|54H must be taken into account when computing 
correlation functions. Thus, the density and field operators become: 



p{x) 
&(x) 



Po(x) + -d x 9(x) 



+ 00 



E 



Po(x) + -d x 6(x) 



m— — oo 
1/2 +oc 



m— — oo 
-■ 1/2 +oo 



p (x) + -d x 9(x) 

7T 



J o " d X po(x') ^2mi9( X ) 
^raiix J* d X ' po(x') ^2mi6{ X ) e ~i4>{x) 

— — o 

e (2m+l)i7T J o dx p (x) e (2m+l)i§{ x ) e ~i4>(x) 



(56) 
(57) 
(58) 



m— — oo 



where 0(x) — Q(x) — irpox. 

The above treatment is essentially correct provided that the external potential represents weak perturbation, which 
has been quantified by requiring that |u ex t(x)| <C /i, for < x < L, and that the external potential is a slowly varying 



function. In terms of its Fourier its transform the latter means that u ex t(9) ~ for \q\ 



> n- 1 



p/hv s . Whereas 



the last condition does certainly hold for a shallow harmonic trapping potential, its effect cannot be considered weak, 
especially near the ends of the atomic cloud where typically u ext (a;) ~ p. Therefore, away from the center of the trap 
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this potential cannot be treated within linear response theory and the above results do not apply. However, one can 
always redo the harmonic-fluid approach by considering fluctuations around a smooth density profile po(x), which 
can be obtained, e.g. from the equation of state in the local-density approximation (LDA) [l£|. Thus, it is not hard 
to see that the effective Hamiltonian becomes: 



h 

2^ 



dx 



vj(x)(d x (j)(x)) 2 + v N (x) (d x @{x) - np (x)Y 



where in the LDA vj(x) = %Trpo(x)/M (i.e. the local Fermi velocity) and 

v N (x) 



1 



On 
~d~ P 



(59) 



(60) 



p=p (x) 



The equation of state p, = p{p) is obtained from the Bethe-ansatz solution [32ll33| . Thus, diagonalization of (|59|l can 
only be accomplished numerically in the general case, and will not be pursued here. 



III. CORRELATION FUNCTIONS 
A. Periodic Boundary Conditions (PBC's) 

We begin by considering the static correlation functions at T = (see Sect. IIII Ul for results at finite temperature) 
for particles in a ring, i.e. periodic boundary conditions. The time-dependent correlation functions can be also 




FIG. 6: Oscillating part (i.e. m — 1 term in jSP ) of the density correlation function as a function of distance, x (normalized 
to the system's size), for several values of K (the different curves are shifted for clarity). As K decreases the system becomes 
more "density stiff", which means that density correlations between distant points are enhanced. However, for large K 
density correlations decay very rapidly with distance. This strong suppression of density fluctuations is characteristic of Bose- 
Einstein condensates in higher dimensions. The dashed line represents the power-law behavior expected for the envelope in the 
thermodynamic limit L — > oo. 
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obtained, but will not be discussed in this paper (except for the dynamic density response function, which is discussed 
in Sect II V C)l . Here we only notice that it suffices to perform a careful analytical continuation to real time t —> it in 
the expressions given in appendices |0 an d El Time-dependent correlations can be important for future experiments 
such like the two-photon Raman out-coupling experiment recently proposed by Luxat and Griffin |18| . 

When computing correlation functions within the harmonic-fluid approach, two strategies are possible. One can 
directly work with the mode expansions given in sections III CI and III Dl using the properties of the exponentials of 
linear combinations of the b(q) and tf(q) operators. Another possibility is to use the more sophisticated techniques 
of Conformal Field Theory (CFT). The latter is our choice here and it is explained in the appendices lei and iDl 
Therefore, the reader unfamiliar with these methods should consult the appendices for full details. In this subsection 
we present the results for the correlation functions for both fcrmions and bosons in a ring. 

Let us first consider the density correlation function, which does not depend on statistics. We first notice that the 
structure of the density operator, Eq. JSj, implies that the density correlation function is a series of harmonics of the 
Fermi momentum pp — irpo. The origin of this has been already discussed in Sect. ITLAl In our calculations, we have 
kept only the leading (i.e. the slowest decaying) term of each harmonic. Thus, 



(p(x)p(0)) = <^{d x Q(x)d x e(0))+p 2 e 2 ™ npoX (A2mAx)A-2 m (0)) pb c 

\ m— — oo 

2m 2 K 



aJi K 



i 



p d(x\L) 



m>0 



1 



p d(x\L) 



cos (2irmpox) 



(61) 



This result is valid in the scaling limit, i.e. for \x\ 3> a, where a is the short-distance cut-off introduced in previous 
sections. The vertex operators A m . n (x,T) are defined in Appendix ICl The coefficients a m are non-universal, in other 
words, they depend on the microscopic details of the model and in general cannot be fixed by the harmonic- fluid 
approach. The function d{x\L) = L\ sm(nx/ L)\/tt is called cord function: it measures the length of a cord between 
two points separated by an arc x in a ring of circumference L. 




FIG. 7: Phase correlations in a ring as a function of distance, x (normalized to the system's size), for several values of the 
dimensionless parameter K. Note that as K decreases the system becomes less "phase stiff", i.e. phase correlations decay 
faster. For large values of K phase correlations decay very slowly, leading to a larger degree of phase coherence between distant 
points. The dashed lines correspond to the power-law behavior occurring in the thermodynamic limit L — > +oo. The final 
upturn of the curves is due to the periodicity of the correlation functions in a ring. 
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We next take up the boson and fermion one-particle density matrices. Upon using Eqs. © and l|10(l . we obtain: 



+00 



(f f (x)*(0)) = po e 2 ™ a; (A 2m ,_ 1 ( a ;)A_ 2m , +1 (0))p bc 



m— — 00 



Pa 



p d(x\L) 



+OO 



p Q d(x\L) 



2m 2 K 



cos {2irmpox) 



(62) 



<*£.(aO* F (0)) =Po E e 2 " (m+ ^ )poX <^2 m+ i,-i(^)A_ 2m _ li+1 (0)) 



pbc 



m=— 00 
+00 

m=0 



1 



p e?(a;|£) 



sin 



27r(m + -)po-^ 



(63) 



Also in these expressions the dimensionless coefficients b m and f m are non-universal. The results of Ref. |30j are 
recovered in the thermodynamic limit L — > +00, which effectively amounts to performing the replacement d(x\L) — > |a;| 
in the above expressions [O . In the Tonks limit if = 1 and the leading term of (|62|l agrees with the exact asymptotic 
results obtained in Refs.|43,|43i 



(*+(a!)*(0))p bc = po65 



Tonks 



1 



-,1/2 



Pod(^l^) 



(64) 



2- 1 / 3 V^A" 6 ~ 0.5214 (A ~ 1.2824271 is Glaisher's constant). 



where &J onks 

In Fig. we have plotted the m = 1 term of (p(x)p(0)) to illustrate how the density correlations are enhanced as 
K is decreased. Thus for K = 1 (free fermions or the Tonks limit of delta-interacting bosons) the distant points are 
more correlated in density than for large K (very attractive fermions or bosons with weakly repulsive interactions). 
This means that the oscillatory terms become more important as K decreases, whereas for large K they can be safely 
neglected, as it is done in the Bogoliubov approximation. On the other hand, as shown in Fig. distant points 
become less phase correlated for small K as the result of density and phase being conjugated fields. The system is 
then said to be less "phase stiff" . However, for larger K phase correlations are enhanced, and only in this sense one 
can speak of "Bose-Einstein condensate" even though, strictly speaking, there is not such a thing as a condensate 
(i.e. long range-order) in one dimension (for instance (e*^^^) p bc = 0). At most, all that exists in the thermodynamic 
limit is an algebraic decay of phase correlations, which is often called quasi long-range order (QLRO). 



B. Open Boundary Conditions (OBC's) 



In this subsection we consider the ground state correlations for particles in a box. The expressions for the correlation 
functions in this case are somewhat more complicated than for periodic boundary conditions. The reason is that in a 
box translational invariance is lost and two-point correlation functions depend on both arguments separately and not 
only on their difference. Furthermore, the ground state expectation value of some operators becomes non-trivial, that 
is, there are non-trivial one-point correlation functions. This is the case of the density operator (see below). However, 
the operator e 1 ^^ continues to have zero-expectation value, implying the absence of any continuous symmetry- 
breaking even though the system is finite and bounded. 

We first compute the ground state expectation value of the density: 



+00 



(p(x))= Po e 2 ' mw (VoW)obc-Po l+E 



m>0 



1 



p d(2x\2L) 



K 



cos(2mirpox + S n 



(65) 



where the coefficients c m and S m are model dependent [Hj]. The above expression is valid in the scaling limit, which 
for the above expression means that min{x, L — x} 3> a, i.e. sufficiently far from the boundaries. It is noticeable that 
(p(x)) exhibits Friedel oscillations, independently of the statistics of the constituent particles. This is not surprising 
in view of the many common features exhibited by interacting bosons and fermions in one-dimension, and whose 
origin was already discussed in Sect. ^ In the thermodynamic limit, the replacement d(2x\2L) — * \2x\ shows that 
the different oscillating terms in the expression for (p(x)} decay from the boundary as power laws with increasingly 
large exponents: m 2 K for m = 1, 2, . . . indicating that only the leading two terms are important. 
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Next we consider the more complicated two-point correlation functions. We begin with the density-density corre- 
lation function, 



1 + °° 

( P (x)p(x')) = {—(d x e(x)d x ,@(x')) obc + z 2 ™ po{mx ~ mx) (A 2m Ax)A-2 m > fl {x')) ohc 



pI i 



K 

2^ 



m,rn' — — oo 
i 2 



pod(x — x f \2L) 



+ 00 

E 



in. in'— x 



d(x + x'|2L) 
z'|2L) 



if 
" 2^2 

2mm' K 



p d(x + x'\2L) 



O 2inpo {rnx—Tn'x') 



[pod(2x\2L)} m2K [p d(2x / \2L)} m ' 2K 



(66) 



Using the results of Appendix iDl one can also obtain the one-particle density matrices of bosons and fermions. The 
corresponding expressions read: 



+oo 



(*t(x)*(x')) = Po E e 2 ^( m - m ' a; ')(A 2mi _ 1 (x)A_ 2 „ iM ( a ;')>obc 



m,m' — — oo 
-1 



Po V d ( 2a; l 2i ) d ( 2a; 'l 2i ) 



d(x + x \2L) 



2K -(-OO 



d(s-^'|2i) 



m,m'=— oo 
p 2i-Kpo (mx — m' x ) 



[p d(2x|2L)] m2 ^[p d(2x'|2L)] m ' 2K 



+ 00 



{tf F {x)^ F {x')) = Po e 2 ""°[( m+ ^ a; -( m ' + ^^(A 2m+1 ,_ 1 (x)A_ 2m ,_ 14 ( a ;'))ob 



in. in — — oo 



(67) 




(x-x')/L (x-x')/L 

FIG. 8: Phase correlations for a system of bosons in a box. In the left plot, x' has been fixed to the center of the box, while 
phase correlations are sample as x is shifted towards the right end. In the right plot, x' is close to the left end and as x increases 
we move towards the right end. As remarked in the case of a ring, the system becomes less phase stiff as K is decreased. This 
means that in this case the boundary conditions have a large effect on the phase correlations in the Tonks limit (K — 1) than 
in the weakly interacting limit (K S> 1). 
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Po 



P Q 1 y /d(2x\2L)d(2x'\2L) 



d(x + x'\2L)d(x-x'\2L) 



+OG 

£ 

m,m / =— oo 



v — i(m-\-m' its gn(x — x' ) /2 

Jm,7n' G 



d(x + x'\2L) 



d{x - x'\2L) 



2(m+i)(m' + i)K 



J2np a [(m+|)a;-(m'+i)a:'] 



[p d(2x\2L)} 



p d(2x'\2L)} 



(m' + i) 2 K • 



(68) 



The coefficients a 



and c m .m' are model-dependent complex numbers. In contrast with the case of periodic 
boundary conditions, where these prefactors can shown to be real, this cannot be established in the present case |70| . 
The above expressions are accurate in the scaling limit, which for two-point correlation functions means that \x — x'\ 3> 
a as well as mm{x,x',L — x,L — x'} a. It is also worth mentioning that the bulk behavior is recovered for 
|x — a;' | <C min{a;, x',L — x,L~ x'}, i.e. mostly near the center of the box (see Fig.[S|). In this limit, the less oscillating 
terms are those with m — to' , which exhibit, in the thermodynamic limit, the same algebraic decay as the correlations 
for the ring geometry. It is also worth pointing out that in the Tonks limit (i.e. for K — 1), the leading term of 
Eq. i|67|) agrees with the asymptotic result obtained in Ref. l43l 



(*t(z)tf(z')>obc 



. .Tonks 
Po b 00 



p^ 1 v / d(2x\2L)d(2x'\2L) 
d(x + x'\2L)d(x - x'\2L) 



(69) 



where b^ nks = 2- 1 / 3 0re A~ & ~ 0.5214 (A ~ 1.2824271 is Glaisher's constant). Note that this is the same prefactor 
as the one found for PBC's, which is not suprising given that both correlation functions have the same bulk limit. 

It is interesting to point out that for particles in a box there are two types of exponents that characterize the decay 
of correlation functions. Consider for instance the density correlation function, Eq. (|66|l . Its bulk behavior (i.e. near 
the center of the box) is governed by the exponent 2m 2 K for a harmonic that oscillates as e ± 2m7r PoO-a: )_ However, 
the asymptotic behavior when one of the coordinates is taken near the boundary is governed by a different exponent. 
Thus if we set x' a and x ^> x' , a term that oscillates as e ± 2im7r Pox m corre l a tion function falls off with the 
exponent equal to m 2 K, i.e. the same exponent that we have already encountered for the Friedel oscillations of the 
density. For the boson density matrix the corresponding exponents are m 2 K + 3/4if whereas for the fermions they are 
(to + 1/2) 2 K + i/AK. It is important to stress that the presence of the boundary breaks the (Lorentz) invariance of 
-ffeff making space and time directions non-equivalent. To see this, consider a dynamic correlation function such like 
the Green's function (^(x, t)^(x', 0)), for bosons, or (^p(x, t)^(x', 0)), for fermions (t > 0). Using the expressions 
of Appendix [Dj and performing the analytical continuation t —> it after setting x = x' i=s a, one obtains that the 
leading term decays as t~ x l K for t a/v s (for both bosons and fermions). This in contrast with the behavior of the 
same functions in the bulk, where space and time are equivalent (i.e . Lorentz invariance is restored) and one expects 
a behavior like t~ x l 2K for bosons and t _ ( 1 / K+Ar )/ 2 for fermions. [75| 

To illustrate both finite-size and confinement effects we show in Fig. [8] the behavior of leading term in Eq. I|67|l . 
which corresponds to the phase fluctuations {e~ l ^ x ^ e 1 ^^) h c (i--e. m = ml = 0). As mentioned above, the correlation 
function is now a function of x and x' separately. Thus, in Fig. [S] we have considered two situations. In the plot on 
the left, we take x' to be the center of the box. Thus it can be seen that the power-law (bulk) behavior is a good 
approximation for x around L/2, whereas the phase correlations deviate from it as x approaches the right end of the 
box (i.e. for x — x' > 0.3L. These deviations become more important as K is decreased and the system approaches 
the Tonks limit (K = 1). This is because, as remarked in the previous subsection, the system becomes less phase stiff 
as K decreases. Thus, for smaller K the boundary conditions have a larger effect in decreasing the phase correlations 
near the ends of the box. However, for larger K the system exhibits a larger phase stiffness, and phase correlations 
decay very slowly, even near the boundaries. Furthermore, the deviations from the power-law behavior are smaller. 
However, these deviations become much larger with distance if now x' — 0.1L, i.e. when x' is fixed near the left end 
of the box. Then, the bulk power-law is no longer a good approximation and we can observe deviations from it for 
a decade: from x — x' ~ 0.1L to x — x' ~ L. The reasons for this deviations is that initially the correlations decay 
according to the bulk power-law, but they rapidly cross-over to a decay governed by the exponent 3/4K (recall that 
to = to' = 0). Near the right boundary, the correlation function must vanish. This is because the pinning of the 
density field Q(x) at x — L makes fluctuate the phase at that point wildly. 



C. Correlations at finite temperature 

As described in Appendix IU1 the same methods that have allowed us to obtain the correlation functions of a finite 
ring at zero temperature also allow to obtain expressions for the correlation functions of an infinite system at finite 
temperature (see below). 
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At finite temperature one finds (see Appendix Q that correlations decay exponentially with distance (as opposed 
to the algebraic decay found at T = for L — > +00). Therefore, in the expressions below, we shall keep only the 
leading terms, which effectively means that we cut-off the series of harmonics at to = (for the boson density matrix, 
and at \m\ = 1 for the density correlation function and the fermion density matrix. Thus the following expressions 
are obtained: 



<P(Z)P(0)>T 

(tft(a:)¥(0)) 
(*£,0»0*f(O)>t 



tt/Lt 



PoC 



po sinh (ttx/Lt) 
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tt/Lt 



B cos{2ttpqx) 



iv /Li 



po sinh (ttx/Lt) 



2 A' 



po sinh (ttx/Lt) 
tt/Lt 



+ . 



po sinh (ttx/Lt) 



2 ' 2K 



sin \TTpox\ + • • 



(70) 

(71) 
(72) 



the dimensionless coefficients A, B, and C are model dependent. In Scct. lIVBl we shall discuss how to fix the prefactor 
for the bosonic density matrix, Eq. J7TJ. 

From the previous expressions one can see that the behavior of the correlations "crosses over" from algebraic 
decay for a <C \x\ < Lt to exponential ~ exp [— \x\/L c (T)] for \x\ 3> L C (T). The characteristic decay length L C (T) 
depends on the correlation function. Thus, for the non-oscillating part of (p(x)p(0))t, L c (T) ~ Lt- However, for 
the oscillating part L C (T) ~ L p (T) — Tiv s /ttKT — Uvn/ttT. For the phase correlations, however, the characteristic 
distance is of the order of L^T) = Tiv s K/ttT = Hvj/ttT = h 2 p /MT, the latter expression being valid only in systems 
with Galilean invariance. One can use the previous expressions also in finite-size systems provided the temperature 
is high enough. The value of T for which these characteristic lengths become of the order of the system size, L, 
give us a rough estimate of the temperature scale below which finite-size and confinement effects are important. 
Since, in general, Lt L^ ^= L p , this temperature scale depends on the correlation function. Thus, for the phase 
correlations, L^(T) rs L for ~ Tivj/ttL = h 2 po/ML, the latter expression applying to systems with Galilean 
invariance. However, for the density correlations, Lt/tt ~ L or L p (T) w L, whichever gives the highest value of T. 
Hence, T p = vaax.{hv s /TrL, Tiv s /ttKL}. 



IV. DELTA-INTERACTING BOSONS AS A LUTTINGER LIQUID 

A. Extracting the Luttinger liquid parameters from the Bethe-ansatz solution 

In what follows, we shall illustrate many of the concepts introduced above by considering a model of bosons 
interacting with a zero-range potential |32t l33| . This is relevant to cold atomic vapors, where atoms mainly interact 
through the s-wave scattering channel, and their interaction is parametrized by the s-wave scattering length, a s . 
Olshanii |(| has shown that the interaction between atoms confined in ID wave guide with transverse harmonic 
confinement is well described by a zero-range interaction potential v(x — x') = g 8(x — x'), where 

„ _ ^ a.//x (73) 



M£ ± 1 - C(a s /£ ± )' 



where a s is the three-dimensional scattering length, C = 1.4603 . . ., and £j_ — ^fi/Muj± is the transverse oscillator 
length. There is a single dimensionless parameter, 7 = Mg/h 2 po, which characterizes the different regimes (some 
authors prefer to use the ID gas parameter po|oid|, which is related to 7 by 7 = 2/ po\a\D\) ■ The model thus defined 
is exactly solvable. Its solution for periodic boundary conditions was found by Lieb and Liniger [32I |3j| , whereas 
the case with open boundary conditions was solved by Gaudin . Lieb and Liniger were also able to compute the 
ground state energy per particle as well as the chemical potential and the sound velocity by numerically solving a set 
of integral equations. One can thus obtain the phase and density stiffness, vj and vn, by writing integral equations 
which allow one to obtain these parameters. These equations were written down by Haldane j5^|. However, we shall 
not follow this path here but instead, we shall adopt a different approach [5^, |5(| . We first numerically find the 
solution of the Bethe-ansatz equations EH , 

No /k-k\ 

kiL = 2ttI, + a - 2 V tan" 1 3 - . (74) 
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FIG. 9: Density and phase stiffness as a function of 7 = Mg/h 2 po. The dots are the result of the numerical calculation of vj 
from the Bethe-ansatz equations. The dashed line corresponds to the Bogoliubov result for vm, Eq. 1761 . The dotted-dashed 
line is the large-7 asymptotic result, Eq. 1771 . 



When deriving these equations we have assumed that the boundary conditions are not periodic but twisted 
as in i|32|) . where a is the twist in the phase. In the ground state the set of integers {ii} j=1 Ng — 
{-{N - l)/2, -(N - 3)/2, .. ., (N - 3)/2, (N - l)/2} jH EM- The ground state energy for N a particles' and a 
phase twist a can be computed from the solution to the above equations, 



E (No,a)=Y,TW ^ 



2M 

and using Eq. (|31|l and l|37|) and vj can be numerically computed [xij. The results of these calculations are 
shown in figures 151 and 1101 Fig. |^1 shows the behavior of the phase (vj) and density stiffness (vn) as a function of 
the dimensionless parameter 7 = Mg/h 2 po. We have numerically obtained vj to show explicitly that, because of the 
Galilean invariance of the model, vj = vp. 

As already pointed out by Lieb and Liniger [32I l33l | , the Bogoliubov approximation for the ground state energy 
yields a very accurate expression of the sound velocity for 7 < 10. Hence, 



V N = V F 

However, for 7 > 1 one can extract an asymptotic expression for vn from either the expressions for \i obtained for 
large 7 by Lieb and Liniger ^ij^ or from a strong coupling expansion of energy [24ll3l| . which yields: 

v N =v F (l- ® J (77) 

In Fig. EH we have plotted the parameters v s and K as function of 7. They also show good agreement with the 
asymptotic results at small 7, 

, s = ^(l_^Y /2 , ( 78 ) 

/ r-.\ -I/ 2 
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TABLE I: Translation table 7 to K for some values of K used in this paper. 



7 


K 


+00 


1 


3.517 


2 


0.7126 


4 



and large 7, 

v s = v F (l - , (80) 
JT-(l + i). (81) 

From Fig. llOl one can see that as 7 varies from zero to infinity, K varies from infinity to 1. Thus, K = 1 corresponds to 
the Tonks limit, which is also the value of K for non-interacting fermions. Finally, in table[I]we provide a "translation 
table" between K and 7 for the values of K used in plots of this article. 



B. Remarks on cut-offs and prefactors 



One of the disadvantages of the harmonic-fluid approach is that it cannot provide explicit expressions for the 
prefactors of the correlation functions. As obtained within this approach, they are cut-off dependent quantities 
and their values depend on the particular regularization scheme used in the calculation This is a signature of non- 
universality, i.e. the prefactors depend on the microscopic details of the model. Thus, they must be fixed by other 
means, either by comparing with numerical results or by adopting a physically sensible regularization scheme (this is 
possible in certain models) and using rather sophisticated form- factor methods (see e.g. 58] and references therein). 
Sometimes expressions for the prefactors can be obtained in certain limits by working directly with the microscopic 
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FIG. 10: Luttinger- liquid parameters vs. 7 = Mg/h 2 pa. The dashed lines are the small 7 approximations obtained from 
Bogoliubov theory whereas the dotted-dashed lines correspond to the asymptotic expressions for large 7. 
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model. This is the case for the Tonks limit of delta-interacting bosons |42j, [59j or for a number of lattice models, such 
like the quantum-Ising model and or the XY model |^| . 

Before jumping into the discussion of the prefactors for the model of delta-interacting bosons, let us pause for a 
moment to consider a question that will be relevant for the fixing of the prefactors, namely the dependence of the 
wave-number cut-off q c = p/Tiv s on 7. As mentioned in previous sections, for delta-interacting bosons the different 
regimes are characterized by a single dimensionless parameter, 7 = Mg/Ti 2 po, where g is the strength of the interaction. 
For weakly interacting bosons (i.e. 7 — > 0) q c — > = 7 1 / 2 /?o, where £ = h/y/M p is the healing length. On the 
other hand, for strongly interacting bosons (7 — -> +00), q c — > npo/2. The full dependence can be obtained from the 
numerical solution of the Bethe-ansatz equations (|?4*|l ■ The result for q c has been plotted in Fig. ^] 

As we have mentioned at the beginning of this section, within the harmonic-fluid approach the prefactors depend 
of the cut-off, a < q~ x (see Appendix Q. In the weakly interacting regime, Popov |2(j, |48| was able to obtain the 
prefactor of the one-body density matrix . In our notation, his result reads: 

/ 2-cy\ \/T7 27r 

9l ( X ) = (^(xMQ))=p (-^j , (82) 

for 7 « 1. In the previous expression C = 0.57721 ... is Euler's constant and £ the healing length, £ = p^ 1 ^ 1 / 2 . 
Taking into account that the numerical factor e 2 ~ c /4 = 1.037 ~ 1, we notice that the above result can be obtained 
within the harmonic-fluid approach by setting the short-distance cut-off a = £. Indeed, this makes a lot of sense since, 
as discussed above, in the weakly interacting limit one expects that a = q^ 1 = £. Things become more interesting 
when one expresses both the exponent and the prefactor in terms of the Luttinger-liquid parameter K, so that one 
can make contact with the way gi{x) was written in previous sections, 

(^(xMO))=Po f° P ^ K \ (83) 

where we have used that, in the weakly interacting limit, K = tx^ 1 ! 2 ^> 1. The dimensionless prefactor (neglecting 
the difference e 2 - c /4 - 1 ~ 0.037) is 

A Popav (K)= (— J . (84) 
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TABLE II: Numerical prefactor of the one-body density matrix as obtained from Eq. 184H 
(DMC) results. 



exact and diffusion Monte Carlo 



7 = Mg/h 2 p 



+00 
2 x 10 3 
6.6667 

2.0 
0.0667 
2 x 10 -3 



K 



K from the exponent (Exact/DMC) 



Ap opov (K) 



A 



Exact/DMC 



1 

1.002 
1.584 
2.523 
12.425 
70.50 



1 

1.0 

1.582 
2.525 
12.469 
70.42 



0.5642 
0.565 
0.806 
0.957 
1.06 
1.02 



0.5214 
0.530 
0.760 
0.951 
1.06 
1.02 




In view of this result, one may be very tempted to use this expression for any value of K relevant to this model 
(i.e. 1 < K < +00). The result, when compared with the diffusion Monte Carlo data reported in Ref. |2l| and exact 
results "0, EH H3 1 i s vei T surprising (see Table ITljl . It is found that Eq. ll%4*|) is accurate within less than 10% over 
the whole range of K values, becoming essentially exact for K ^> 1, as expected. This fact is quite remarkable taking 
into account that from K « 10 to K = 1 the prefactor changes by a almost a factor of two. 



C. Tips to compare with the experiments in cold gases 



In the previous sections we have equipped ourselves with the tools to characterize the behavior of one-dimensional 
cold-atom systems. In this characterization the Luttinger-liquid parameters v s and K play an prominent role. The 
latter is especially important because it enters the exponents of the correlation functions, and governs their algebraic 
decay at zero temperature. Furthermore, the value of K also tells us to which kind of instabilities the system can 
become unstable when perturbed. A well-known example |ltl l3fj| | is the quantum-phase transition to a Mott-insulator 
in the presence of a commensurate periodic potential, which takes place (strictly speaking at T = and L — > 00) 
for arbitrarily weak potentials at K < 2. At finite temperature, the behavior of the system is characterized by the 
correlation lengths Lt, L^iT) and L p (T), which are combinations of the parameters v s ,K and the temperature, T. 
As we are interested in the bulk properties of ID systems, we shall focus in this section on the implications of the 
correlation functions previously derived for the experimental observations. 
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We begin our discussion by considering the finite temperature momentum distribution. This quantity can be 
accessed by performing Bragg scattering measurements at large momentum transfer [Tsl l67l l68[ . At temperatures 
T > T = Ti 2 p a /ML, one can safely neglect finite-size and boundary effects. In this temperature range the effects of 
the inhomo gen eity in harmonically trapped systems can be treated using the local density approximation as described 
in Refs.l67ll68l For a uniform system, the momentum distribution is defined as the Fourier transform of g\(x,T), i.e. 

p + CO 

n(p,T)~ / dxe~ ipx g x (x,T) (85) 



After introducing Eq. (|71|) into this expression and approximating the prefactor by Ap opov (K) = {K/ir) 1 t 2K , we 
obtain the following expression: 



T (1-1/ AK + ipL^(T)/2K) 



(86) 



where L^T) — h 2 po/MT = poA 2 -, (At = h/s/MT is the de Broglie thermal wave-length) is the phase correlation 
length. This expression for n(p, T) is different from the lorentzian form commonly used throughout the literature. 
The lorentzian results from taking the Fourier transform of 



5l(X ' T) " P °(d^) 



1/2K 



-\x\/2L$(T) 



(87) 



which is the asymptotic form of g\(x,T) for |x| 3> Lt- Hence 

'2\ Tk 



In Figs. [T21 andUlwe compare both forms, Eq.® and JBEI>; n(p, T)/n(0, T) and n'(p,T)/n'(0,T) are plotted 
against pL<f,(T). Our calculations assume the scaling limit, which in these plots means that \p\ < q c . For \p\ 3> q c , 
Olshanii and Dunjko have shown |19j that n(p, T) ~ p~ 4 . From figures lT2*lll3l and 1 141 one can draw the conclusion that 
the lorentzian is a good approximation for large K, i.e. for weakly interacting bosons. The width of the lorentzian 
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FIG. 14: Momentum distribution at finite temperature for a system with K — 4. The inset shows that the asymptotic behavior 
is not as p~ 2 , but rather as p 1 ^ 2K ~ 1 = p~ 7//8 . 



is proportional to L^,(T) oc po/T, which thus provides a measure of this ratio. However, in the Tonks limit the 
momentum distribution at T > looks more like a stretched lorentzian. The reason is that, in the large momentum 
limit, n(p,T) as given by Eq. I|87|) . behaves as a power-law: p _ ( 1_1 / 2K ) ^ p -1 / 2 , for K = 1. Experimentally, 
extracting the Luttinger-liquid parameter from n(p,T) may be hard because the true n(p,T) may not display this 
algebraic behavior in a sufficiently wide range of p (set aside the complications that averaging g\ using the local density 
approximation may introduce). However, there may be some chance to observe the more qualitative behavior shown by 
these plots: as the system's parameters are tuned (e.g. by using a Feshbach or the confinement-induced resonance Q) 
towards the Tonks limit, the momentum distribution should become more "stretched" and less Lorentzian-like. 

A more quantitative estimate of the Luttinger-liquid parameters near the center of the trap can be obtained from 
the dynamic density response function, 

X (q,u>) = [ dxdte iut -^ X (x,t), (89) 



X (x,t) = -jrO(t)(\p{x,t),p(0,0)]) T , (90) 

whose imaginary part is measured in low-momentum Bragg scattering experiments |44| . In the low momentum (= 
long wave-length) limit, we can replace p(x,t) by IL(x,t) — d x 6(x, t)/ir. Thus, using the equations of motion for 
II(x, i), one can obtain the result (rj — > + ): 

X(<l,u) = — . 27T, (91) 

nh (id + %f]Y — Lj z (q) 

where ui(q) = v s \q\ for |g| <§; q c . Hence, 

2 

Im X (q, to) = [S (w - u(q)) - 5 (« + w(g))] . (92) 

Since vj — v s K, the weight becomes Kq/2. Thus by measuring the dispersion oj{q) one can obtain v s , whereas from 
the weight of lmx(q, w), K can be measured. Inhomogeneity effects can be treated also in this case within the local 
density approximation, as discussed in Ref . I44L In this respect, we notice that Imx(g,w) has the same structure as 
the result from Bogoliubov approximation in d > 1. The only difference is that it holds, as well as the expression for 
n(p, T) derived above, for all values of the interaction parameter 7, and not only in the weakly interacting limit. 
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V. JASTROW-BIJL WAVE FUNCTIONS AND LUTTINGER LIQUIDS 



In Sect.[n]it has been shown that the low-energy spectrum of a Luttinger liquid is completely exhausted by phonons 
with linear dispersion. This phonons are nothing but collective oscillations of the density and the phase. In view of 
this fact, one may wonder what kind of correlations are built into the ground state wave function by these collective 
excitations. The question has been considered in different contexts by several authors 0, ^(|. Here we restrict 
ourselves to the case of particles in a box (OBC's). Taking inspiration from the work of Reatto and Chester we 
write the ground state wave function for Nq bosons in a box as follows: 



$0(^1, . . . , x N ) = $ c (xi, . . . , x N )$ (xi, . . . , x N ). 



(93) 



where 3> c (xi, . . . , Xn) describes the correlations between the particles in the ground state whereas &o(xi, . . . , xn) 
describes their independent motion, which for No bosons in a box is given by 



No /nx ., 
$0(2:1, . . .,x No ) = J^sin (-^-J 



(94) 



that is, all the Nq lie in the lowest-energy orbital sin(7ra;/L). As for $ c (xi, . . . ,xn ) we will obtain an asymptotic 
expression valid for \xi — Xj\ 3> a, for i,j = 1, . .., N . In this limit the correlations between particles are dominated 
by the zero-point motion of the collective excitations. 

In order to obtain the asymptotic behavior of <I> C , we consider the ground-state wave function of the low-energy 
Hamiltonian (in this section A denotes the operator of the corresponding classical variable A) , 



H c ff 



hv. 



dx 



^ii 2 (x) + * (dj(xj) 



(95) 



For OBC's, it is convenient to use the following expansions: 



n(x) 

4>(x) 



n„ 



q>0 
9>0 



(96) 
(97) 



where q — rnr/L, with n a positive integer and IIo = N — Nq . The canonical commutation relations now read 

n o ,0o = i, 

tin. d>n' — iS„_„'. 



(99) 



We shall work in first quantization, which means that the wave function $ c is regarded as functional of H(x). From 
the second of the above commutation rules 



. d 



(100) 



Therefore, when expressed in terms of Fourier modes the Hamiltonian takes the form 



hnv s ^ 



2K 



q>0 



qK Yx 



frirvs \ -\ 
2K ^ 

q>0 



n- 



qK 



dlP q 



(101) 



i.e. it represents a collection of decoupled harmonic oscillators, and therefore its ground state is just a product of 
gaussians: 



exp 



2K ^ q 

q>0 y 



(102) 
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When expressed in terms of II(ie), this function becomes: 



$ c = cxpji^ dxj dx'U(x)K(x,x')Il(x')\ , 



(103) 



where 



/C(x, x') 



E 



-«</ 



2n 



LK ^ q 

q>0 y 



■ cos(<7x) cos(qx') 



K 



In 



4 sin 



n{x + x 
2L 



n(x — x') 
2L 



(104) 



Notice that K(x,x') depends on the short-distance cut-off a and therefore, it cannot correctly describe short-distance 
correlations where x ~ x' . The second expression above is thus only asymptotically correct. 

As we discussed in Sect.|n] the configurations of H(x) describe long wave-length density fluctuations, i.e. the "slow" 
part of the density. Thus, we are quite tempted to make the replacement H(x) — ■> p(x) — ^<=i $( x i ~ x j)- However, 
in doing this we must be careful enough to remove the terms that involve the function K(x,x / ) evaluated at x = x' 
since K (x, a;') cannot describe short-distance correlations (In Ref. ^3 we failed to notice this, and as a result the wave 
function obtained there is not completely correct. See below for the correct expression). If we do so, the following 
function is obtained in terms of particle coordinates: 



$c(xi, . 



,XN ) 



n 



Sill 



n(xi + Xj) 



2L 



sm 



ir(xi - xj) 



2L 



2N (N - 
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TTX, 
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cos 
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l/K 



(105) 



Hence, the complete ground-state wave function reads 

®o(xi, . . -,x No ) = AfJJ cos (^j 1 ) - cos ( ^ l 



1<J 



l/K 



No 

II< 

i=l 



V L 



(106) 



where J\f is the normalization constant. 

There are several reasons to believe that the wave function (|106fl is, at least asymptotically, correct. First, in the 
non-interacting limit K — > +oo (i.e. 7 — * for delta-interacting bosons), $ c — ► 1, and one recovers the independent- 
particle ground state $0 = $o- On the other hand, in the Tonks limit K — 1, $0 is the exact ground state of a system 
of hard-core bosons in a box ,43|. Furthermore, if one repeats the above calculation for bosons in a ring, the result is 
the well-know Jastrow-Bijl function: 



$™ s (x, 



,Xn ) 



TTiXi 



L 



It turns out that this is the exact ground state wave function of the Calogero- Sutherland model 
of hard-core bosons with long-range interactions, and whose Hamiltonian reads [7^: 
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It is also known that this model also has a solution under harmonic confinement 



i.e. for 



(107) 

which is a model 
(108) 



1 M^^, 2 n 2 K^jK- 1 -!) 

cs ~ 2M ^ dxj + 2 + M ^ ( Xi - Xjf ' 

i=l 1 i=l i<j v 



The exact ground-state wave function has a similar structure to Eq. (|106|) : 
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There is also a good chance that the latter wave function may be asymptotically correct for other models of interacting 
bosons in a harmonic trap. 
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VI. CONCLUSIONS 



In this paper we have addressed the properties of one-dimensional systems of cold atoms using the harmonic-fluid 
approach (also know as "bosonization"). Besides reviewing the method in pedagogical detail, we have argued that 
it allows to treat boson and fermion systems both in strongly and weakly interacting limits. We have also shown 
how concepts and results obtained using the Bogoliubov-Popov approach 0, US and its modifications 0, |2(| can be 
naturally recovered using bosonization, which in our opinion is much simpler conceptually When combined with the 
conformal field theory methods explained in the appendices, it becomes a very powerful tool to obtain the functional 
forms of the correlation functions for various geometries (ring, box with Dirichlct BC's), and also finite-temperature 
correlations. The method has some limitations, however, as it cannot provide explicit expressions for the prefactors 
of the correlation functions, which turn out to be model-dependent. Moreover, for strongly interacting systems it 
becomes difficult, without further input, to relate the phenomenological parameters vj and vn to the microscopic 
parameters of the model at hand. Nevertheless, for a relevant model of bosonic cold atoms, we have shown that 
one can successfully extract those parameters from the exact (Bethe-ansatz) solution. Furthermore, the prefactor of 
the one-body density matrix could be fixed with less than 10% error by expressing a result obtained in the weakly 
interacting limit [48j in terms of K = y/vj/vjy. These results should allow for a quantitative comparison with the 
experiment. In this respect, we have discussed how to extract the Luttinger-liquid parameters K and v s = ^vj^vj 
(i.e. the sound velocity) from Bragg-scattering measurements of the density response function and finite-temperature 
momentum distribution. For the latter quantity we have argued that, as a system is tuned into the Tonks regime, 
the finite-T momentum distribution should become more stretched, as compared to the lorentzian form exhibited by 
n(p, T) in the weakly interacting limit. Looking forward to experiments where these predictions can be tested, we 
hope that the present work will foster the use of the harmonic-fluid approach in the study of ID cold-atom systems. 
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which states that 11(2:) and the phase 4>(x) are canonically conjugated in a low energy subspace. We now proceed to 
give a heuristic derivation of it (a different argumentation can be found in Appendix [5] m terms of path integrals) . 
Our starting point will be the following commutation relation between the density and momentum density operators: 
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APPENDIX A: COMMUTATION RELATIONS OF U(x) AND 4>(x) 



In the main text we have often used the commutation relation: 



[U(x),<j)(x')] = iS(x-x'), 



(Al) 



\fi(x),j p (x')} = iKd x >S(x - x')p(x'). 
This result can be most easily derived by working in first quantization, where 



(A2) 



N 




(A3) 




(A4) 



i=i 
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In the previous expressions Xi and pi = —ihd Xi stand for the position and momentum operators, respectively. Next, 
we shall consider long wave-length fluctuations of the density and the current. Thus we set p(x) ~ po + II(ie) and, 
upon linearizing, 

j p (x) = * [&(x)d x *(x) - d x tf(x)*(x)] « n Po d x cj>(x), (A5) 

we arrive at 

[H(x),d x/ (f>(x')] = id x ,5{x - x 1 ). (A6) 
By integrating this expression over x' , it reduces to 

[U(x),(p{x')]=iS(x-x') + C. (A7) 
The integration constant C must be zero since from this commutation relation it must follow that 

[N,e-*+W] =e-^ x \ (A8) 

which follows from the fact that the field operator ty'(x) adds particles to the system. Thus we have provided a 
justification for Eq. IjAlf) . Alternatively, one could have written Eq. I|A6() as 

[d x Q(x), d x i(f>{x')} = md x i8{x — x') = —ind x S(x — x'), (A9) 

which upon integration over x becomes: 

[0(a), d x 4{x')\ = -iwS(x - x') + C (A10) 

To show that the constant C = one has to work a little bit more in this case. For open boundary conditions one 
has to work out this result for < x,x' < L by using the mode expansions, Eqs. (|47I48|I . and then formally taking 
the limit a — > + . For periodic boundary conditions, however, this is required for the expression: 

[j,e ]=i (Aii) 

to hold. Thus, after setting C = and defining the canonical momentum Tl^lx) = d x (f)(x)/7t, we can write the above 
commutation relation as: 

[IL <t ,{x),e(x')] = id(x-x'). (A12) 

This is just a different representation of the duality of the fields 0(x) and 4>(x), which provide two complementary 
descriptions of the low-energy physics. Finally, we urge those readers unhappy with the rather non-rigorous treatment 
of operators in this appendix and Sect.[n]to consult appendix iBl to be reassured of the correctness of the results. 

APPENDIX B: PATH INTEGRAL FORMULATION 

Our goal in this appendix is to present a somewhat different, but at the same time complementary derivation of 
some aspects of the harmonic-fluid approach. To this purpose, we will employ the path integral formalism. We also 
compute the low-temperature limit of the partition function to show that the spectrum described by this formulation 
has the same structure and degeneracies as the one obtained from the Hamiltonian given in Eq. (|28|l of Sect. Ill CI 

As explained in e.g. Ref. l69l the partition function of a bosonic system in the grand canonical ensemble, Z can be 
written as a coherent-state path integral: 



[dip*dip] e ~ s[rM . (Bl) 
The functional S[ip* , ip] is the (euclidean) action, and for the the Hamiltonian in Eq. I|ll|l has the following form 



1 



S^*,^] = / — I dx hip*(x,T)d T tp(x,T) - nip* (x, t)iP(x,t) + —\d x ip(x,r)\ 



2m 

+ 77 / dx' v(x - x)i/;*(x,t)iP*(x',t)iI)(x',t)iI)(x,t) . (B2) 
2 Jo 



1 ' L 
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Using "polar coordinates", ip(x,r) = \fp{x, T)e^^' r ) and ip*(x,r) — p(x, r)e ^( X > T ) , where p(x,r) and 4>(x,t) are 
real functions (i.e. not operators), the action becomes: 



fit* 



o h Jo 



I dx 



ihp(x, r)d T (j){x, r) + —p(x, r) (d x </>(x, t)) 2 + -d T p(x, r) 
2m 2 



d x p(x,r) 
8m \ p{x,t) 



H p(x, t) 



1 



dx' p(x, t)i>(x — x')p{x' , r) 



(B3) 



Next we proceed to give a more explicit meaning to the coarse-graining procedure used in Sect.[Hl As explained there, 
we first split: 



p(x,r) = p<(x,t) +p>(x,t), 
(f)(x,r) = (t)<(x,T) +cp > (x,r), 



(B4) 
(B5) 



where p>(x, r) and </>>(x, r) describe the fast modes, i.e. those with momenta higher than a -1 = min{i? _1 , q c }, and 
frequencies higher than oj c ~ /x. The fields p<(x, r), ^<(x, r) describe the slow modes. To give a more mathematically 
precise definition, we can define the slow part of a given field, h{x, r) = p(x, r), <p(x, r) as 

h<(x, t) = J dx'dr' f(x -x',t- T ')h{x', r'), (B6) 

where f(x,r) is a slowly varying function over distances ~ a and imaginary times ~ h/u> c . The low-temperature 
description of the system is obtained by coarse-graining the action, i.e. by integrating out the fast modes. Thus we 
define the effective-low energy action by 

e- s ^ e ^ - J [dp>#>] e -^>.*>.e.*], (B7) 

where, just as we did in Sect, [n] we have parametrized the slow modes by <f> = 0<(x, r) and 0(x,t), such that 
p<(x, t) = 9a;0(x, r)/7r. In the general case, performing the functional integral involved in Eq. (|B7|) it is not feasible. 
However, on physical grounds and from the structure of the microscopic action, Eq. (|B3|) . one can guess the following 
form [7U: 



dr I dx 
'o 



-d x Q{x, t) d T 4>(x, t) + {d x (j)(x, r)) 2 

7T Z7T 



(9 x 9(x,r) - 7r/9 ) 2 
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(B8) 



For weakly interacting systems one can perform the coarse-graining perturbatively , and to the lowest order this 
amounts to keeping only the quadratic terms in the gradients of 8 and 0. Thus, for instance, for delta-interacting 
bosons one has vj = v s K — Tmp^/M = vp, as discussed in Sect.|n] and = vJK = g/lvK + 0(g 2 ). For this model, 
one can also show, using the effective fermionic Hamiltonian reported in Ref. |23, 0] j that the action has the form 
(|B8ll . with vj — vf and vn = vp (l — 87 _1 + 0(7~ 2 )). Such a derivation will be given elsewhere j6|| (the result for 
vn can be also obtained from the expressions given by Lieb and Liniger for the chemical potential at large 7). 

Note that the imaginary term (i.e. the Berry phase) ihd x O d T (f)/ir indicates that d x Q/-K (and hence H(x,t)) is 
canonically conjugated to <p, just as in the original action the Berry phase reflected that %j> and ip* are canonically 
conjugated fields (see e.g. Ref. 169^ . 

For bosons, the path integral must be performed over configurations that obey ip{x, t + Tij3) = ip(x, r) besides the 
periodic boundary conditions, ip(x + L,t) = ip(x, r). Hence, from the expressions for the field operators in terms of 
O and 4>, these fields must obey the following boundary conditions: 



Q{x + mL,T + nTil3) = Q(x, r) + muN + mrP, 
<p{x + mL.T + nhf3) = 4>{x, t) + ran J + nirQ, 



(B9) 
(BIO) 



where m,n are arbitrary integers, and (— I)' 7 = (—1)^ = +1. 

After obtaining the low-energy effective action, S e g, we notice that it describes a quadratic theory. Therefore, we 
have two choices: to integrate out 9 or to integrate out tfi. This yields two apparently different representations of the 
theory. If we integrate out <p and introduce 6 = — Trpox we obtain 



1 



5cff [ 



2-kK 



dr I dx 
Jo 



- {d T ef + v s {d x ef 



(Bll) 
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Had we integrated out instead, 



S e $[4>] 



K 



hp ph 
dr I dx 
Jo 



— (d T (j)) 2 + v s {d x <t>) 2 



Sb- 



(B12) 



where Sb — ipo J dx [<p(x,h/3) — 4>(x,0)j. This is a "surface" term that can be in principle dropped. However, we 
shall retain it in the computation of the partition function that follows because we are working with a finite-size system. 
Notice that both representations, Eqs. i|Bll|) and (|B12|) . are dual to each other, in the sense that one can be obtained 
from the other by means of the replacements 6 4> and K <-> K . This property has important consequences, in 
particular the correlation functions of one field alone can be obtained from the corresponding correlation function of 
the other field alone by means of the replacement K — > K . Another important feature of these representations 
is that they show that the role played by K is that of an effective "temperature" or Planck's constant: When K 
is large, the density field, 9, fluctuates wildly whereas the phase field, (j>, behaves almost classically. For instance, 
this is what happens in a system of weakly interacting bosons (a "quasi-condensate" at finite T, which resembles 
more and more a BEC as the temperature is lowered). On the other hand, when K is small it is the phase that 
fluctuates more violently whereas the density behaves more classically. This corresponds to a system whose behavior 
approaches crystallization. However, in both cases no symmetry breaking can take place in one-dimension, as both 
phase transitions (BEC and crystallization) break continuous symmetries. At most, as it is emphasized in the main 
text, one gets quasi-long range order. 

We finally compute the partition function using the representation of the action in terms of the phase field, Eq. (|B 1 2|) . 
Upon writing (j)(x, r) = 4>o(t) + irxJ/L + -KTQ/(Tif3) + <p(x, r), where (— 1) J = (— 1)^ = 1, and <Po(t + fi(3) = 4>o{ T ) is 
the spatially homogeneous part of 4>{x, r) , whereas 







u>,q=£0 



i{qx — WT) 



<j>(q,uj) 



(B13) 



with q = 2nm/L and u> — 2im/ {%[}), and m ^ 0, n integers. At this point it must be stressed that these Fourier 
expansions must be cut-off at a momentum q ~ a . Effectively, this can achieved by restricting ourselves to 
temperatures much lower than Hv s /a < /i. 

Using the above decomposition, the partition function can be written as the product: 



Z — ZqZqZjZ . 



(B14) 



where 



Z 

Zq 
Zj 

z 



2tt 

#o(Q) 



KL 



exp • 



2ttv, 



h6 



dr 



s JO 
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inNoQ 



[cu 2 + (v s q) 2 ] 



(B15) 
(B16) 
(B17) 

(B18) 



<t£0,uj 



Notice that Zq is formally equal to the imaginary-time propagator of a particle of mass Mq = KL/(jrv s ) moving 
on a line, subjected to the boundary condition that after a time t = h/3 it should return to the starting point: 
4>(h(3) = (f>o (0). From this observation and the form of the imaginary-time propagator in one dimension: G(r) = 
y/M /2TrT exp{-M (^(r) - 0o(O)) 2 /2r}, we obtain that 



dM0)G(hp) 



<2KL 

hf3v s 



(B19) 



The integral over the initial position < <^>o(0) < 2tt stems from the fact that the partition function is a trace. 
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Next consider the other terms of the product. Using Poisson summation formula we can rewrite Zq as follows (we 
set Q = 2q): 



X>p{-^f } c2 ""'"- = J^JZ iz exp {~W^k + 2 " (JV " N ° )z } <™> 
£ ^{-^wt^ -"<■?}• < B21 > 



* a +°° 
2KX 

JV=— oo 

Thus we can define 

+00 



Z N ,j = Z Z N Zj= £ exp{-/3^^(7V_7V ) 2 } £ exp A ■ (B22) 

N=—oo ' even ,7 ^ ' 

The fact that vY runs over all integers follows from the requirement that Q is even, which is related to the bosonic 
statistics of the particles, as it has been discussed above. 

Finally, we recognize in Z the partition function of a system of non-interacting bosons 69] with linear dispersion 
uj(q) = Hv s \q\ > 0, 



-1 



Z = \{ M -e-^'Wj . (B23) 
Hence, 

Z= £ ^{-^(N-N A £ expj-^^^ln^-e-^^)" 1 . (B24) 

N=-oo ' even J ' <j^0 

This result is the same that can be obtained from the Hamiltonian in Eq. I|28() . This shows that the path integral and 
Hamiltonian formulations describe the same spectrum with the same degeneracies. 



APPENDIX C: FINITE SIZE/TEMPERATURE CORRELATION FUNCTIONS 

In this appendix we shall compute two important correlations functions, from which any two-point correlation 
function can be derived. Before we do it we will introduce a lot of new technology. These methods are related to a 
symmetry of the phonon model defined by H c g known as conformal invariance, which has deep consequences. This 
is a too vast subject to be covered here, and we refer the interested reader to the literature [5Sl 17(1 l7l| . However, 
in the discussion below we shall try to be as self-contained as possible. The introduction of this technology will be 
lengthy at the beginning, but it really makes life much easier when complicated correlation functions with different 
boundary conditions need to be computed. This means that part of the results of this appendix will be also used in 
the following one. Here we begin by considering a system with periodic boundary conditions. Let us first define the 
vertex operators: 

A m , n (x,T) = e *™e{*,T) e iM{*,T) i ( C1 ) 
One of the correlation functions in which we are interested is the following 

(Am,n(x, t)A— n (x, T)) p bc- (C2) 

The other one is 

{d x e(x,T)d x ,e( x ' y)) phc . (C3) 

We have explicitly written that these correlation functions obey PBC's to emphasize this aspect with respect to other 
correlation functions that will show up below. In the above expressions r stands for the imaginary time (see below). 

Before we proceed any further, it is useful to introduce a new set of fields, 4>l{x) and <pi,(x), implicitly defined by 
the following expressions: 

1 [feW-^Wl, (C4) 



2VK 
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We next employ the mode expansions presented in Sect. Ill CI It is convenient to work with operators that depend 
on the imaginary time r. To find the dependence on t of the fields 0# and 4>Lt one solves the equations of motion for 
the normal modes. For instance, 



^ & (^ T ) = [H e s,b(q, t)] = -hv s \q\b(q,T) 



which yields b(q,r) = b(q)e v *\i\ T t etc. When this is also done for N, J as well as for <f> and 8 , one finds 



2wi ^— \ / zvi 
4>l{x,t) = -j-N L (v s T + ix) + } I — 
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i\X,T) = +(j} R + — N R {V S T - IX) + I ~T 

where q = 2nm/L {in — 0, ±1, ±2, . . .) and the zero modes: 

(pQL = VK(po 
4>or = VK(j) 
N L 
N R 
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,-aq/2 



b{-q)e- q(VsT+lx) + b\-q)e q{v > T+ix) 
b{q)e- q{VsT - 1 ^ + tf{q)e q{VsT - lx) 
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(C6) 

(C7) 
(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 
(C13) 



Note that [N L ,<p 0L ] = [N R ,<t>o R ] = i but [N L ,<j> m ] = [N Rl (j) 0L \ = 0, as follows from [N,cj> ] = [J,0 ] = i. Thus, 
4>r{x,t) and 4>l{x,t) commute, which implies that the operator 



A mn (x t) = e ^( m '~ n )^- L ( :r ' T ) e W m > n )0«O r > T ) 
where /3(m, n) = m^/K /2 + nj2\[K. Furthermore, 

H ^ = ^ f dx [(dMx)) 2 + (d x Mx)f 



(C14) 



(C15) 



Thus, the Hamiltonian splits into two independent parts. The reason for this is that </>i(x,r) and 4>r{x,t) represent 
modes propagating in opposite directions. This can be more clearly seen by making an analytic continuation to real 
time: r — > it. Thus we see that <pL{x,t) = 4>l{x + v s t) and therefore describes the modes propagating to the left 
(hence the subindex L), whereas <f> R (x,t) = <ft R (x — v s t) represents the modes propagating to the right (hence the 
subindcx R). This property is called chirality. 

To make contact with conformal field theory it is useful to introduce complex coordinates w — v s t + ix and 
w = v s t - ix, and to denote z = e 27Tw/L and z = e 27T ™ /L . Thus, 



0l(x,t) = (j> L (z) = ~0ol ~ iN L \nz + V" —= 

m—l v 

+00 ^ 

4>r(x,t) = <j> R (z) — +4>qr + iN R \nz + V] — = 

i V m 

m—l v 

We also introduce the chiral vertex operators: 

Vp(z) — : e il3 ^ L<yZ ' > : — e ~ i f 3 't> aL e +P NL lnz 

Vp{z) -- 



2nm 
2ixm 



+ z m 6 f 
+ z m 6 t 



2nm 



exp 



+ 00 



e i(34> R (z) . _ e if3<pQ Re -!3N R \nz Q ^ 



m—l 



27rm\ 

L K 


exp 


m—l v v 7 


'irm\ 
L ). 


exp 


m—l v x 7 





(C16) 
(C17) 

(C18) 
(C19) 
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where : . . . : means the operators are normal ordered as indicated above. However, the vertex operator A ntm (x, r), as 
given by Eq. (j(J14|) . is not normal ordered. To write it in normal order form one must take into account that: 



e i04> L (z) = a p J /2 



— j V (z), (C20) 

e M R W = a ff/2 f?™y 12 (C21) 



where a is the short-distance cut-off introduced in Sect III Cl However, in what follows we shall not keep track of 
those factors and will simply use the replacements e l ^^ L ^ — > Vp(z) and e*^ R ^ — > Vp(z) in Eq. I|C14|) . As described 
below, these factors can be restored at a latter time by performing a conformal transformation. 

The vertex operators just introduced have an interesting property that makes easier the computation of n— point 
correlation functions. This property is: 

{VpM — V^Zp)) = Hizj-zuf^, (C22) 

j<k 

(V 01 (z 1 )---V 0v {z p )) = H( Z] - Zk p^ (C23) 

j<k 

provided that X^i=i Pj = ((• " ') stands for the expectation over the ground state of H e g, see below). This is some 
times called the neutrality condition. It can be regarded as a consequence that the Hamiltonian, Eq. (|C15|) is invariant 
under the infinitesimal shifts <f> <l — > 4> 'l + e and 4>r —> <j)R + e. Another way of proving this condition is to notice that 

p 

(cxp{ij2 [pim^-nMoL +^K,™ J )fe]}) = (e'^r^V'^!"^") (C24) 
i=i 

vanishes unless ■ mj — Yl j n j — because the ground state is an eigenstate of J and N. At this point it is important 
to remark that in the calculations that follow we assume that in the ground state (J) =0, (N) = Nq (No being integer, 
which implies that we work in the canonical ensemble). This is the case for bosons, but for fermions the selection rule 
(— l)' 7 = — (— 1) N requires N to be odd. However, if N is even, then the ground state is degenerate and (J) = ±1. 
The implications of this degeneracy are discussed below. The neutrality condition as stated in Eq. (|C24(1 implies that 
the only allowed vertex operators in this field theory are those for which p = f3(m, —n) for the left-moving field and 
(3 — f3(m, n) for the right-moving field, m and n being integers. Other choices lead to non-integral exponentials of 8q 
and 4>q, which in general are not physical. 

Next we shall prove Eq. (|C22|I . We shall first do it for two vertex operators (i.e. p = 2) and later explain how to 
extend the proof to p > 2 (we consider left- moving fields only, the proof for right-moving fields is identical). Thus, 

{Vp (zi)V-p (z 2 )} = (e- 113 ^ eP Nt - ln Zl e Woi e~ pNL ln 22 e'^ {zi ] e l ^ L {zi ] e~ ip *L We"** 1 ) . (C25) 

Using the identity e A e B = e^ A ' B ^e B e A and [$ L (zi), $^(z 2 )] = — ln(l - Za/zi), where <J>l(z) = 
Sm=i z ~ m h (—2nm/L) / \fm, one finds that 

(VpizxW-pifr)) = e -/3 2 ln Zle -/5 2 ln(l-, 2 / Z1 ) = ^ _ ^-^ (C26) 

Thus the rule to compute the expectation value of an arbitrary number of vertex operators is to first separate the zero 
modes from and <&\. Next, commute all the terms involving and 4>ol to the left of the terms involving $l and 
Nl- Every time one does so with a pair of operators i and j , one gets factor exp \ftiPj ln(zj — Zj)] = (zj — z^^i . The 
final expression is thus rendered into normal order, that is, it has all the operators Nl to the right of 3>t , </>oi> 
The expectation value of the normal ordered product equals unity provided the neutrality condition is obeyed (and 
the ground state is not degenerate). 
Using the above results, 

(A m ,n{z-L,Zi)A- m - n (Z2,Zl)) = (yp(j n> - n ){zi)V-p( m< - n ){Z2)) (Va( m ,n) (^l)^-/3(m,n) (#2)) 

= (zi- z 2 )- p2 ^ 1 --^ (zi- Z2 )-P 2 (™< n ) (C27) 
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which follows upon using Eq. I|C14|I and equations l|C22IC23|) . However, this result does not yet produce the correct 
correlation function for PBC's by simply replacing z — ► e 2 n{v 3 T+ix)/L an( j ^ _^ g 2n(v s T-ix)/L^ b ecause f ^ ne factors 

thrown away after normal ordering A mn (x, r). The way to recover those factors is to regard the replacement 



z = e 2 ™' L 



z = e 2 ™' L 



(C28) 



as a conformal transformation that maps the infinite complex plane onto an infinite cylinder of circumference L. 
In conformal field theory |70l l7l| , it is known that such transformations lead to multiplicative renormalization of 
the correlation functions of (primary) operators. Thus if Oi{z,z) is a set of (primary) operators with the following 
two-point correlation functions in the coordinates z, z: 



(C29) 



where hi and hi are called conformal dimensions, then under an arbitrary conformal transformation z — z{w) and 
z = z(w), their n— point correlation function transforms as [70| 



P /j \ -hi 

aw 



{O ll (w 1 ,wi) . . .O lp (wp,w p )) = (O il (zi(wi),z 1 (w 1 ))...O ip (Zp(w p ),z p (wp))) J| ( — 



dw \ h% 



(C30) 



It is not hard the see that when the transformation is given by Eq. I)C28|1 . the terms (dw/dz)~ J , etc. give us back 
the factors in Eqs. (|C20IC21|> . Therefore, by applying the above formula to the chiral vertex operators, one obtains 
the following expression: 



(W 2 , W 2 ))pbc = 

tt/L 



( 



2 7TL -l e TT(w 1 +w 2 )/L^I 32 ( m <- n ) / 2nL -l e n(w 1+ w 2 )/L\ P 2 (rn,n) 



sinh (ir(wi — W2)/L) 



ir/L 



sinh(7r(wi — W2)/L) 



- -(- Jl r 

' 4 K 



sinh (n(wi —W2)/L) 



sinh {-k{wi — W2)/L) 



The constants c mn are cut-off and model dependent. Hence, upon setting w\ = ix and u>2 = ix', 



(C31) 



{^m,n (x)-A— m, — n (j£ ))pbc 



l 



pod(x — x'\L) 



^imn-Ksgn(x— x ! ) /2 



(C32) 



where we have introduced the mean density po to get a dimensionless prefactor c m ,„; the cord function d(x\L) = 
L \ sin(7rx/i)|/7r. 

As we have pointed out above, for fermions the ground state is degenerate if (TV) = Nq is even. In this case 
(J) = ±1, which using (|C11(I and (|C11D implies that (Nr) = —{Nl) = \f~K/2(J). Therefore, the expectation value of 
fully normal ordered products of vertex operators no longer equals unity. Instead, for a two-point correlation function 
one is left with the factor 



i)Nl In fi 



w -f3(rn.n)N R ln Jl \ 

)(e y ' ' z ? ) = exp 



(3(m, -n)(N L ) In — - /3(m, n)(N R ) In 2L 

Z2 Z 2 



(C33) 



Upon making the replacements z\ % = e 27r ( VsT1 - 2+lXl - 2 " L and Z\$ — e 27T ^ Vsri - 2 42; i. 2 )/ L , one obtains the following phase 
factor: 



F m , n (xi - x 2 ,n - r 2 ) 



, - t < J) ( mv J ( T l - T 2 ) - in {x i -x 2 ) ) /L 



(C34) 



which after analytical continuation to real time (r — ► it) becomes F m , n (x,t) = e m { J )i mv J' t nx))/L_ Thus, the left- 
moving part of the field operator W F (x), tpr, ~ A+i t —i must multiplied by the factor F+x t -x(x, t) — e - 47r < J )( t, J t + :r )/ L j 
whereas the right-moving part ipu ~ by F_x,-x(x, t) — e m { J ){vjt-x)/L _ When one works in the grand canonical 

emsemble similar phases will appear for both fermions and bosons because (N) — iVo(/i) is not an integer. In both 
cases the phases are of 0(1 jV) and simply disappear as L — > +oo. However, in finite systems they must be taken into 
account if detailed comparison with, e.g. numerics, is being sought. In what follows, however, we shall assume that 
(J) = and Nq is a positive integer (odd in the case of fermions). 
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We next turn to the computation of the second correlation function mentioned above, namely 

(d x 9(xT)d x >9(x' , r')) p bc- We first notice that d x = i (d — d) , where d = (v~ 1 d T — id x )/2 and d — (v~ 1 d T + id x )/2. 
Therefore, we first compute: 

Wz(«iWz(aia)> = -~ 1 —2, (C35) 



(d<pR(zi)d(f> R (z 2 )) 



(zi - z 2 ) 
1 



(C36) 



(z\ - z 2 ) 

where the derivatives are taken with respect to z and z. Using Ea.f IC5)l we write 6 in terms of and tp^. Hence 

1 



K K 

(de(z 1 ,z 1 )d9(z 2 ,z 2 )) = —(d(/)L(zi)d(f> L (z2)) = — -5, 

4 4 [zi - z 2 ) 



(d0(z 1 ,z 1 )B9(z 2 ,z 2 )) 



K 



(d(j> R (zi)dcj) R {z 2 ) 



K 



1 



4 {zx - z 2 y 

Applying the conformal transformation (|C28(I to these correlation functions, one obtains: 

(d6(wi , w 1 )d9{w 2 , w 2 )) p bc 



K 



tt/L 



- K 
(d6(w 1 ,w 1 )d9(w 2 ,w 2 )) phc = - — 



sinh (n(wi — w 2 )/L) 

TT/L 



sinh (jr(wi — w 2 )/L) 



Hence, 



{d x 0(w 1 ,w 1 )d x 0(w2,w 2 )) v toc 



-((d - d)8(w 1 ,w 1 )(d - d)9{w 2 ,w 2 )) phc 
-(dd(w 1 ,w 1 )dd(w 2 ,w 2 )) phc - (d0(w 1 ,w 1 )d9(w 2l w 2 )) phc 



K 

T 



ir/L 



sinh (n(wi — w 2 )/L) 



tt/L 



1 2 



sinh(7r(wi —w 2 )/L) 



In the static case, one sets w\ = ix and w 2 = ix' , and the above expression reduces to: 



{d x 6{x)d x 6{x')) = - 



Kpl 



1 



pod(x — x'\L) 



(C37) 
(C38) 

(C39) 
(C40) 

(C41) 
(C42) 

(C43) 



Finally, we shall describe how to obtain correlation functions at finite temperature in an infinite system. In the 
path integral formulation discussed in Appendix [S] finite temperature means that one has to sum over configurations 
of the bosonic fields satisfying ip(x, r + h/3) = i/j(x, r) (for fermions the boundary conditions are anti-periodic 69]). In 
terms of the complex variables w = v s t + ix and id = v s t + ix, this means that tp(w + Lt, w + Lt) = w), where 
Lt = fiVgfi = fiv s /T is the thermal length. If we compare this with the requirement of PBC's: ip(w + iL, w — iL) = 
ip(w,w), we see that the replacements iL — > Lt, for the left-moving part, and — iL — > Lt, for the right moving 
parts, lead to the sought correlation functions at finite temperature. Otherwise, performing the following conformal 
transformation: 



Z — e 2 ™ w / L T £ _ e -2i-rrw/L T 



on Eqs. 1C27|) and (|C37IC38|I leads to the same results (w = v s t + ix and w = v s t — ix): 



1 (0,0))t = c m ,n 
K { 



tt/Lj 



sin (ttw/Lt) 

tt/Lt 
sin (ttw/Lt) 



(A mjn (w,w)A-. 

(d x 6(w,w)d x 6(0,0)) T = 
Upon setting w = ix and w = —ix, one obtains: 

(A m , n (x)A. m ,- n (0)) T = c m , n j m ™ ssa(x)/2 



tt/L 7 



sin (ttw/Lt) 
tt/Lt 



sin (ttw/Lt) 
sin (ttw/Lt) 



sin (ttw/Lt) 

tt/Lt 



(d x 8(x)d x 0(O)), 



Kpl 



tt/L'_ 



p \ sinh (ttx/Lt) 

2 



2 \_po sinh (ttx/Lt) 

Thus we see that, at finite temperature, the correlation functions fall off exponentially with distance. 



(C44) 

(C45) 
(C46) 

(C47) 
(C48) 
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APPENDIX D: CORRELATION FUNCTIONS WITH OPEN (DIRICHLET) BOUNDARY CONDITIONS. 

In this Appendix we heavily rely on the results of the previous Appendix. We first make the important observation 
that for open boundary conditions the density and phase fields can be expressed in terms of a single chiral boson field, 

<M X ' T ) = 7T7t7 ['^- r ( x ' r ) + ( I } r(~ x i t )] (Dl) 
zv A 



9(x,r) = 9 b + ^[4> r (x,t)-M-x,t)}. (D2) 

The main difference with the chiral (right-moving) field of the previous section is that the field 4>r(x) must be now 
defined for — L < x < +L. Therefore, it will be a function of z = exp [2tt(v s t — ix)/2L]. In terms of this complex 
coordinate, 

+00 

Mz)=<t>OR+iN R \nz+ -7^ [z- m b(^)+z m tf (™)] . (D3) 

m— 1 v 



The zero mode [Nr, 4>or] = i since Nr = (N — No)/yK and 4>qr — yK <po- The reason why only one chiral field is 
need is because with OBC's the wave number q of the excitations can be only positive. Furthermore, the presence of 
boundaries makes impossible the existence of persistent currents, and thus only N can appear in the mode expansions. 
Therefore, only one chiral field operator can be constructed out of (t>o : N and b(q) and b'(q). Such a field, however, 
must be defined for — L < x < L, such that q is quantized as 2nm/2L = mir/L, as required by OBC's. We chose 
the interval < x < L to correspond to the actual system, and below we shall assume that x > always. As in the 
previous Appendix, we can write: 

A m n(x 7") = e im(> ( x > T ) e in <t>( x > T ) = e imS B e i0{m,n)4>R{x,r) e ip{—m,n)4> R {—x,T) (D4) 



where f3(m,n) = {m\JK + n/\/K)/2. Notice as well that since 4>r(—x,t) = 4>r(z*), the fields 4> an d are linear 
combinations of ^>h(z) and 4>r(z*). At this point, it is worth pointing out that in the present context z* 7^ z because 
z and z are to be treated as independent variables, whereas here z* denotes the complex conjugate of z. As we did 
above, we shall replace A m , n (z) by its normal-ordered form: 

Am >n (z,Z*) = V r /3 ( m , n )(z) %(- m , n ){z*)- (D5) 

As a consequence, 

( z P> z p)) - const, x (Vp( — m p ,n p ) 

That is, for OBC's any n-point correlation function of the fields 9 and 4> becomes a 2n-point correlation function of 
the chiral field 4>r. This a particular version of a more general result in boundary conformal field theory obtained by 
Cardy |7l|. In the last expression the neutrality condition still holds. This fact can be again regarded as a consequence 
of the invariance of the Hamiltonian 



#cff = I dx [8 x <t> R (x)] 2 (D7) 



under the infinitesimal shift 4>r — ► (f>R-\-e. It is worth to pause to examine the consequences of the neutrality condition 
for OBC's. The condition requires that 

p 1 p 

J2 [0( TO i> n i) + P(- m i> n i)] = 7T7f Yl n i = °- ( D8 ) 
3=1 3=1 

Notice that now the condition does not involve mi,ra2, . . . ,m p , and therefore does not forbid the existence of non- 
vanishing A m fi(z), but implies that (A nim (x, r)) b c = 0, for any n / 0. Mathematically, this is expected since 9 is 
not affected by the shift 4>r — > 4>r + e or, in other words, it does not contain <f>Q. Therefore, 

_- _ __ m 2 JtT 

(A mt0 (z,Z*)) =C m (^( m ,o)(«)^(-m,0)(«*)) =Cm{z-Z*) * . (D9) 

Next, to obtain correlation functions at finite size, we perform the conformal transformation 

- = e 2™/(2i) = e ir*/L ( D10 ) 
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which maps one half-plane onto a strip of width L. The transformation law for chiral (primary) operators is: 



(Oi^m) . ..O ip (w p )) = (Oi^ziim)) . . . O lp (z p (w p )) J] 

i=i 



dw 

dz 



(Dll) 



where the operators have the two-point correlation function (Oi(z)Oi(O)) = z 2hi . Using the above expression 
(recalling that w = v s t — ix) 



(A mfi (w,w*)) ohc = (A TO>0 (a;)) o bc = c„ 



tt(w+w*)/ L 



/L 



oTTW j L p7TKJ* / L 



4 


1 


— Cm 


_p Q d(2x\2L)_ 



(D12) 



where the mean density po has been introduced to make the pre-factor c m dimensionless. 
Next we compute the two-point correlation function of the vertex operators: 

x(zl — z^P(- m - n )P( m ' ~ n ) (f 2 — 2*^/3(-™,-«)/3( m '>-«) 

After performing the conformal transformation (jDlOfl . it takes the form 



(D13) 



(A mn {wx)A- m :- n {w2))ohc = C n m ,r. 



p Q d(wi - w{\2L) 



4 4K 



1 



Pod(w 2 - w^\2L) 



4 4K 



[pad{wi - wl\2L)] 



(mm'K-^ r )/2 



[podim-w^L)}^" 1 '^^/ 2 



sinh[7r(u;5; - w 2 )/2L] 



sinh[7r(ii;i — w^/^L] 



n(m~m') /4 



sinh[7r(ii;^ — W2)/2L] 
sinh[7r(uii — W2)/2L] 



n(m+m')/4 



(D14) 



In the main text, we set W\ — —ix and u>2 — —ix' so that 



{A mn (jJC^A_ rri/ / n (x ))obc ^m,m' 



c" , 


l 


4 4K 


1 




p d(2x\2L)_ 




[p d(x'\2L)\ 



4 4K 



[p d(x + x'\2L)} {mm ' K -^ )/2 
[p d(x-x'\2L)} (mm ' K+ ^ )/2 



jin(m-\-m )sgn(x — x )/4 



In the above expressions we have assumed that < x = — Im W\,x' = — Im u>2 < L. 
Finally we need compute the following correlation function 

{d x 6{x)d x ,6{x')) ohc . 



d x 6(x, t) = — - [d x <j) R (x, t) + d x <f> R (-x, r)] = — - \d(f) R {w) + d(j) R (w*)] 
2 Li 



To this end, we first notice that Eq. (|D2|) implies 



Using 



we obtain 
K 



(d(j> R (z 1 )dc/) R {z2)) = 



1 



(2l - Z2? 



(2*) s 



K 



(z 1 -z 2 ) 2 {z{-z*) {zl-z 2 ) (zi-z*) 



The conformal transformation in Eq. (|D10(1 along with the transformation property (|D 1 1|) . then yield 



(d x 9(w 1 ,wl)d X '8(w 2 ,W2))obc = -£< 



7T/2L 



sinh(7r(it;i — w 2 )/2L) 

77 /2L 



it/2L 



sinh(7r(u;:j; — u^^L) 



sinh(7r(w^ — iu 2 )/2£) 



it/2L 



1 2 



sinh(7r(w^ - w%)/2L) 



(D15) 

(D16) 
(D17) 
(D18) 
(D19) 

(D20) 
(D21) 
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Upon setting w\ 



-ix and wi = ix', we finally obtain 

{d x e{x)d x ,e{x')) ohc = I 



d(x-x'\2L) 



d{x + x'\2L)_ 



(D22) 
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